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FOREWORD 


This  report  contains  the  abstracts  and  viewgraphs  of  the  presentations  at  the 
Fifvfnth  Anniiai  Mechanics  of  Composites  Review  sponsored  by  the  Materials  Laboratory. 
Each  was  prepared  by  its  presenter  and  is  published  here  unedited.  In  addition^  a 

LISTING  OF  BOTH  THE  IN-HOUSE  AND  CONTRACTUAL  ACTIVITIES  OF  EACH  PARTICIPATING  ORGANI¬ 
ZATION  IS  INCLUDED. 


The  Mechanics  of  Oompositfs  Review  is  designed  to  present  programs  covering 
ACTIVITIES  THROUGHOUT  THE  UNITED  STATES  AlR  FoRCEy  NaVY,  AND  NASA.  PROGRAMS  NOT 
COVERED  IN  THE  PRESENT  REVIEW  ARE  CANDIDATES  FOR  PRESENTATION  AT  FUTURE  MECHANICS  OF 
COMPOSITES  REVIEWS.  ThE  PRESENTATIONS  COVER  BOTH  IN-HOUSE  AND  CONTRACT  PROGRAMS  UNDER 
THE  SPONSORSHIP  OF  THE  PARTICIPATING  ORGANIZATIONS. 


Since  this  is  a  review  of  on-going  programs^  much  of  the  information  in  this 

REPORT  HAS  NOT  BEEN  PUBLISHED  AS  YET  AND  IS  SUBJECT  TO  CHANGEj  BUT  TIMELY  DISSEMINATION 
OF  THE  RAPIDLY  EXPANDING  TECHNOLOGY  OF  ADVANCED  COMPOSITES  IS  DEEMED  HIGHLY  DESIRABLE. 

Works  in  the  area  of  mechanics  of  composites  have  long  been  typified  by  disciplined 

APPROACHES.  It  is  HOPED  THAT  SUCH  A  HIGH  STANDARD  OF  RIGOR  IS  REFLECTED  IN  THE  MAJOR- 
ITY>  IF  NOT  ALL^  OF  THE  PRESENTATIONS  IN  THIS  REPORT. 


Feedback  and  open  critique  of  the  presentations  and  the  review  itself  are  most 

WELCOME  AS  SUGGESTIONS  AND  RECOMMENDATIONS  FROM  ALL  PARTICIPANTS  WILL  BE  CONSIDERED 
THE  /PLANNING  OF  FUTURE  REVIEWS. 


GEORGE  E.  HUSMAN>  Chief 
Nonmetallic  Materials  Division 


Materials  Laboratory 


IN 


ACKNOWLEDGEMENT 


We  express  our  appreciation  to  the  authors  for  their  contributions  and  to  the 

POINTS  OF  CONTACT  WITHIN  THE  ORGANIZATIONS  FOR  THEIR  EFFORTS  IN  SUPPLYING  THE  PROGRAM 
LISTINGS. 


ix 


DETECTION  OF  FAILURE  PROGRESSION  IN  CROSS-PLY  GRAPHITE/ EPOXY 
DURING  FATIGUE  LOADING  THROUGH  ACOUSTIC  EMISSION 


Jona than  Awerbuch ,  William  F,  Eckles  and  Eliezer  Katz 


Department  of  Mechanical  Engineering  and  Mechanics 
Drexe 1  Univ  er s i ty 
Philadelphia,  Pennsylv  ania  19104 


ABSTRACT 


Monitoring  acoustic  emission  (AE)  du 
practical  procedure  for  detecting  fatigue 
destructive  tool  is  particularly  attractive 
acquisition  of  data  In  real-time,  Its  po 
progression  and  accumulation,  for  anticlpa 
different  failure  mechanisms  and  determining 
ty  to  non-v isual  damage . 

Problems  remain  to  be  solved,  however: 
data  obtained,  the  appropriate  test  methodol 
AE  results  and  the  actual  deformation  charac 
tion  of  the  various  failure  mechanisms  an 
emission  generated  by  damage  and  that  genera 

In  this  study,  AE  is  monitored  in  a  var 
of  different  stacking  sequences  and  conta 
Three  loading  sequences  are  applied:  mon 
quasi-static  loading-unloading  while  Increme 
fatigue  loading,  all  In  uniaxial  tension, 
of  nondestructive  (visual,  X-radiography ,  ac 
and  destructive  (laminate  deplylng,  photomlc 
techniques.  AE  is  also  monitored  with 
different  failure  mechanisms  can  be  Isolated 


ring  fatigue  loading  appears  to  offer  a 
damage  and  damage  growth.  This  non- 
because  of  the  simplicity  in  its  use,  the 
tentall  for  monitoring  damage  initiation, 
ting  failure  sites,  for  Identifying  the 
damage  criticality,  and  for  its  sensitlvl- 

proper  interpretation  of  the  voluminous 
ogy  to  be  employed,  the  correlation  between 
teristics  and  state  of  damage,  identlfica- 
d  processes,  and  the  distinction  between 
ted  by  friction . 

lety  of  cross-ply  graphite/epoxy  laminates 
Inlng  different  ratios  of  ply  thickness, 
otonlc  quasi-static  loading  to  failure, 
ntally  increasing  the  load  to  failure,  and 
The  AE  results  are  compared  with  a  variety 
ousto-ultrasonics ,  and  frequencey  response) 
rography,  and  scanning  electron  microscopy) 
specially  designed  lay-ups  in  which  the 


Results  indicate  that  stacking  sequence,  and  to  a  lesser  degree  the  ratio  of  ply 
thicknesses,  strongly  affects  the  event  intensities.  Emphasis  is  placed  on  determining 
the  correlation  between  the  AE  results  and  actual  failure  processes  for  different 
laminates.  During  both  quasi-static  and  fatigue  loading,  a  significant  amount  of 
emission  is  generated  by  friction,  in  some  cases  exceeding  that  generated  by  new 
damage.  During  fatigue  loading,  the  friction  generated  emission  can  be  discriminated 
by  the  load  level  at  which  it  occurs  and  through  proper  correlation  among  the  AE  source 
intensity  variables.  Results  indicate  that  matrix  failure  and  fiber  breakage  result  in 
middle  and  high  range  AE  source  intensities.  Due  to  the  rapidity  with  which  matrix 
cracks  and  delamination  progress,  distinction  of  these  processes  is  more  difficult  and 
investigation  continues  in  this  direction.  Frequency  response  of  the  subject  laminate 
is  shown  to  have  potential  for  determining  the  state  of  damage.  Initiation  of 
transverse  cracks,  damage  accumulation,  and  the  state  of  transverse  crack  saturation 
are  clearly  distinguished  by  this  method  In  cross-ply  laminates. 


CONCLUSIONS 

1.  Stacking  sequence  and  ply  thickness  have  significant  effect  on  damage  initiation 
and  accumulation  and  on  the  failure  process  as  detected  through  acoustic  emission. 

2.  A  significant  amount  of  emission  is  generated  by  friction  among  newly  created 
fracture  surfaces  during  both  quasi-static  and  fatigue  loading. 

3.  The  friction  generated  emission  can  be  distinguished  from  emission  generated  by  new 
damage  through  proper  correlation  among  various  AE  source  Intensities. 

4.  Preliminary  results  indicate  that  AE  source  Intensities  of  friction  generated 
emission  depend  on  the  type  of  failure. 

5.  Monitoring  AE  during  fatigue  loading  can  Indicate  the  cycle  number  at  which  damage 
initiates  and  progresses  and  the  type  of  damage. 


1 


6.  The  frequency  response  of  the  subject  laminate  strongly  depends  on  the 
damage.  Amplitude  reduction  and  shift  In  frequency  can  be  qua  1 1 la t Iv a ly 
with  state  of  damage. 


DETECTION  OF  FAILURE  PROGRESSION  IN  CROSS-PLY 
GRAPHITE/EPOXY  DURING  FATIGUE  LOADING 
THROUGH  ACOUSTIC  EMISSION 


Jonathan  Awerbuch,  William  F.  Eckles  and  Eliezer  Katz 


Department  of  Mechanical  Engineering  and  Mechanics 
Drexel  University 
Philadelphia,  Pennsylvania  19104 


Program  sponsored  by  AFWAL/FDL,  F336 15-84-3204. 
George  P.  Sendeckyj  of  AFWAL/FIBEC  is  the  program 
monitor. 


OBJECTIVE 

DETERMINE  THE  VALIDITY  OF  THE  .ACOUSTIC  EMISSION  TECHNIQUE 
IN  DETECTING  DAMAGE  INITIATION  AND  ACCUMULATION  AND  IN 
IDENTIFYING  THE  FAILWIE  MECHANISMS  AND  PROCESSES  IN 
GRAPHITE/EPOXY  LAMINATES  DURING  FATIGUE  LOADING 


APPROACH 


•  IDENTIFY  THE  APPROPRIATE  TESTING,  DATA  REDUCTION 
AND  DATA  ANALYSIS  METHODOLOGIES 

•  PERFORM  EXPERIMENTAL  PROGRAM  WITH  A  VARIETY 
OF  CROSS-PLY  LAMINATES  AND  DETERMINE  THE  EFFECT 
OF  STACKING  SEQUENCE  AND  PLY  THICKNESS  ON  ACOUSTIC 
EMISSION  RESULTS 

•  DETERMINE  THE  ACOUSTIC  EMISSION  SOURCE  INTENSITIES  WHICH 
BEST  IDENTIFY  MATRIX  CRACKING,  DELAMINATION,  FIBER 
BREAKAGE  AND  FRICTION 

•  CORRELATE  THE  ACOUSTIC  EMISSION  RESULTS  WITH  OTHER 
NONDESTRUCTIVE  (Visual,  X-Radiography,  Acousto-Ultrasonics, 
and  Frequency  Response)  and  DESTRUCTIVE  (Laminate  Deply, 
Photomicrography,  and  SEM)  TECHNIQUES 


State  of 
cor rela  ted 


2 


Figure  3 
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Figure  5. 


Distribution  of  AE  source  Intensities  for  events  accumulated 
during  the  Initial  part  of  fatigue  loading  (R  =  0.1)  of 
graphl te/epoxy  [02/902/0]g  laminate.  Significant  emission  is 
generated  by  friction  among  newly  created  fracture  surfaces. 
High  source  intensities  occur  only  at  the  upper  load  range. 


evehts  events  events  events 


SPEC.  NO.  UN/4  2c/W  »  0.304  N  =  0  -  250  Cycles 

LAY  -  UP:  [0],.  o,  =  122.9  MPa  f  »  0,1  Hz  E  -  1,002  EVENTS 
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EVENTS  EVENTS  EVENTS  EVENTS 


Figure  8. 


events  events  events  events 


SPEC.  NO 


2c/W  =  0.  304 
(7  -  122.9  MPa 


N  =  210  -  250  Cycles 
f  -  0.1  Hz 


AMPLITUDE  (dB) 


AMPLITUDE  (dB) 


ENERGY  COUNTS 


ENERGY  COUNTS 


DURATION  [ysec] 


DURATION  [ysec] 


LOAD  RANGE  =  (95-100%)a, 
E  =  17  EVENTS  ^ 


LOAD  RANGE  =  (10-60%) a, 
E  =  131  EVENTS 


Figure  9.  Distribution  of  AE  source  intensities  for  events  accumulated 
during  cycle  210  through  250  of  the  same  specimen  shown  in 
Figure  7  for  a)  events  generated  at  the  upper  part  of  the 
load  range;  b)  events  generated  at  the  lower  part  of  the 
load  range.  Identification  of  damage  through  AE  may  require 
analysis  of  several  source  intensities  and  load  ranges. 


SPEC.  NO.  5C1/2 


for  graphite /epoxy  [02/902/0]g  laminate ;  a)  unloaded  specimen;  b)  specimen 
under  75%  of  maximum  previous  load  cycle.  Amplitude  at  frequencies  of  150 
and  470  KHz  reduces  dramatically  with  damage  initiation.  In  the  frequency 
range  of  0.8  -  1.0  MHz  amplitude  is  constant  but  a  shift  in  frequency  occurs 
during  damage  accumulation.  Following  dama ge  saturation,  frequency  is 
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ABSTRACT 

The  composite  laminates  now  being  used  in  aircraft  and  aerospace  structups  are  sub¬ 
jected  to  a  variety  of  impact  threats,  such  as  dropped  tools  during  fabrication  and 
handling,  or  hard  objects  such  as  runway  stones,  ice  balls,  etc.  Research  results 
during  the  past  two  decades  indicate  that  such  impact  damage  is  frequently  non-visual 
but  nevertheless  results  in  significant  performance  degradation.  This  sensitivity  of 
composites  to  non-vlsual  Internal  damage  Is  of  major  concern,  primarily  when  these 
materials  are  used  in  primary  structures. 

Among  the  conventional  nondestructive  techniques  for  detecting  and  locating  non-yisual 
Impact  damage  are  ultrasonic  C-scan  and  X-radiography .  Monitoring  acoustic  emission 
(AE)  is  a  promising  NDT  procedure  in  that  it  can  provide  information  about  damage  pro¬ 
gression  and  damage  criticality  in  real-time.  Consequently,  this  study  focuses  on  the 
applicability  of  the  AE  technique  to  impact-damaged  composites  subjected  to  external 
loading . 

This  research  constitutes  a  comprehensive  program  to  provide  a  detailed  characteriz¬ 
ation  of  Impact  damage  In  graphl te /epoxy  laminates  subjected  to  low  velocity  (20-500 
m/sec)  Impact  and  to  develop  the  proper  procedures  foE  monitoring  acoustic  emission  In 
Impact-damaged  laminates.  A  variety  of  experimental  and  analytical  techniques  are  em¬ 
ployed  In  order  to  determine  Impact  damage  criticality,  with  primary  emphasis  being 
placed  on  the  applicability  of  the  acoustic  emission  technique  to  detect  and  locate 
non-vlsual  Impact  damage,  monitor  Its  progression  under  external  loading,  and  deter- 
mine  its  criticality.  Four  major  areas  are  addressed  in  this  program:  1.  nonde¬ 

structive  detection  and  evaluation  with  emphasis  on  the  acoustic  emission  technique; 
2.  destructive  examination  for  detailed  mapping  of  the  impact-damaged  region;  J. 
failure,  fracture,  degradation  in  performance  and  their  modeling;  and  4.  determina¬ 
tion  of  damage  criticality  to  identify  the  major  causes  of  degradation  in  properties 
and  the  NDT  detection  threshold. 

The  acoustic  emission  results  Indicate  that  artificially  Induced  damage  (e.g.  circular 
holes  and  silts)  could  be  easily  detected  and  located  at  relatively  early  stages  of 
loading.  Similar  results  were  obtained  with  laminates  containing  complete  perforation 
due  to  Impact.  However,  when  non-vlsual  Impact  damage  occurs,  with  significant  delam- 
Inatlon,  matrix  cracking  and  splitting  (as  detected  through  stereo  X-radlographs ,  ult¬ 
rasonic  C-and  F-scans,  photomicrographs  and  deplled  laminates),  detection  and  location 
of  existing  damage  during  proof  loading  Is  more  difficult,  mainly  because  of  the  sign- 
ificant  emission  generated  from  throughout  the  specimen  length  during  the  rapid  pro 
gresslon  and  accumulation  of  the  matrix  dominated  failures.  It  has  been  determined, 
through  the  analysis  of  the  AE  event  Intensities,  that  during  quasl-statlc  loading 
significant  emission  Is  generated  by  friction  between  Impact-damaged  fracture  surfaces 
In  contact  with  each  other.  The  more  severe  the  damge  Is,  the  larger  the  amount  of 
friction  generated  emission  and  the  lower  the  stress  level  at  which  It  Initiates. 
Only  at  higher  loads  Is  actual  damage  progression  detected  through  the  AE.  Conse¬ 
quently,  an  experimental  procedure  and  AE  data  analysis  methodology  have  been  develop¬ 
ed  to  detect  and  locate  non-vlsual  Impact  damage  during  cyclic  loading.  The  results 
Indicate  that,  based  on  the  friction  generated  emission,  such  damage  could  be  eas^y 
detected  and  located  with  a  few  load  cycles  and  for  dynamic  stress  levels  less  than 
30%  of  static  ultimate  stress.  With  an  Increasing  number  of  cycles,  damage  Initiation 
and  progression  could  be  Identified  as  well. 

other  NDT  techniques  and  examination  procedures  developed  and  applied  In  this  research 
program  Include  the  acous to-u 1 trason Ic  (AD)  technique  and  the  frequency  response  of 
the  Impact-damaged  laminate.  The  results  Indicate  that  both  techniques  can  easily 
detect  the  existence  of  damage  and  Its  severity  and  they  qualitatively  correlate  well 
with  the  strength  degradation  of  the  subject  laminate.  From  stereo  X-radlographs  the 
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depth  of  the  delaminated  interfaces  and  the  extent  of  splitting  could  be  approximately 
Identified,  and  these  results  also  correlated  well  with  the  results  obtained  with  the 
ultrasonic  F-Span  and  deplylng  techniques.  Other  issues  addressed  in  this  study 
Include  the  effects  of  artificially  Induced  damage  (e.g.  delamination,  broken  fibers 
and  notches),  oblique  impact,  laminate  thickness  and  configuration,  strength  degrada¬ 
tion,  etc.,  most  of  which  are  discussed  in  the  formal  presentation. 

CONCLUSIONS 


1.  Impact  damage  is  largely  non-visual,  nevertheless  nondestructive  and  destructive 
examinations  reveal  significant  Internal  damage  in  the  form  of  de lamina tion , 
intraply  cracks,  matrix  splitting  in  all  damaged  plies,  and  fiber  breakage.  e- 
gree  of  stiffness  and  strength  degradation  depend  on  laminate  configuration. 

2.  During  quasi-static  loading  in  tension,  the  detection  and  location  of  non-visual 
Impact  damage  through  acoustic  emission  is  difficult.  Significant  emission  is 
generated  by  friction  among  the  impact-damaged  fracture  surfaces. 

3.  However,  based  on  this  friction  generated  emission,  during  cyclic  loading  impact 
damage  could  be  easily  detected  and  located  with  a  few  load  cycles  and  at  rela¬ 
tively  low  dynamic  stress  levels.  Damage  progression  could  be  easily  tracked  and 
failure  mechanisms  identified. 

4.  Damage  severity  could  be  easily  determined  through  the  acous to-ul tra sonic  tech¬ 
nique.  All  event  intensities  strongly  depend  on  the  level  of  damage  and  rapidly 
attenuate  with  increasing  severity. 

5.  Similar  results  are  obtained  using  the  frequency  response  of  the  subject  lami¬ 
nate,  which  was  found  to  be  highly  sensitive  to  impact  damage.  It  seems  that  the 
frequency  response  depends  on  the  details  of  the  Impact— damaged  region. 

6.  Good  qualitative  correlation  could  be  established  between  the  various  destructive 
and  nondestructive  examination  procedures. 
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OBJECTIVES 


•  DETERMINE  THE  APPLICABILITY  OF  THE  ACOUSTIC  EMISSION 
TECHNIQUE  TO  DETECT  AND  LOCATE  IMPACT  DAMAGE  IN 
COMPOSITE  LAMINATES,  TRACK  ITS  PROGRESSION,  IDENTIFY  THE 
FAILURE  MECHANISMS,  AND  EVALUATE  ITS  CRITICALITY 

•  CORRELATE  THE  ACOUSTIC  EMISSION  RESULTS  WITH  THE  MOST 
COMMONLY  USED  NONDESTRUCTIVE  AND  DESTRUCTIVE 
EXAMINATION  TECHNIQUES 

•  CONDUCT  A  DETAILED  CHARACTERIZATION  OF  THE  IMPACT- 
DAMAGED  REGION  AND  IDENTIFY  THE  MAJOR  MATERIAL 
VARIABLES  AFFECTING  ITS  SEVERITY 


APPROACH 


1.  DETECTION  TECHNIQUES  II. 

EXAMINATION  TECHNIQUES 

(NON-DESTRUCTIVE) 

(DESTRUCTIVE) 

1.  VISUAL 

t. 

PHOTOMICROGRAPHY 

2.  X-RADIOGRAPHS  (3-D) 

2. 

DEPLYING 

3.  ULTRASONICS  (C  and  F-  SCANS) 

3. 

S.E.M.  (3-D) 

4.  ACOUSTO  -  ULTRASONICS 

5.  FREQUENCY  RESPONSE 

6.  ACOUSTIC  EMISSION 

7.  C.C.T.V. 

III.  FAILURE  AND  FRACTURE 

1.  STIFFNESS  AND  STRENGTH  DEGRADATION 

2.  APPLICATION  OF  FRACTURE  MODELS  APPLIED  TO  COMPOSITES 

3.  FINITE  ELEMENT  SIMULATION  OF  IMPACT  DAMAGE 

IV.  DAMAGE  CRITICALITY 

1.  CORRELATION  WITH  ARTIFICIALLY  INDUCED  DAMAGE 
(Broken  Fibers,  Inserted  delamination.  Notches) 

2.  OBLIQUE  IMPACT 

3.  STACKING  SEQUENCE  EFFECT(S) 

4.  LAMINATE  CONFIGURATION  (EFFECTS) 

5.  LOADING  FUNCTION  EFFECT(S) 

6.  THICKNESS  EFFECT(S) 
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SPEC.  NO,  CrOI-21  SPEC,  NO.  BQI-1  SPEC.  NO.  HQI-13 

V.  =  39.6  m/s  V.  -  135.6  m/s  V,  =  39.3  m/ s 

=  732.6  MPa  -  292.6  MPa  =  191.1  MPa 

E^=  19,863  EVENTS  E  =  2,697  EVENTS  =  0.1  Hz 
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AMPLITUDE  (dB) 


a  =  292.6  MPa 
=  2,406  EVENTS 


a.  SPEC.  NO.  GQI-21 

V.  =  39.6  m/s 
1 


a  =  732.6  MPa 
=  25,650  EVENTS 


86.0  100. 


b.  SPEC.  NO.  BQI-1 
V^  =  135.6  m/s 

so  ■  I  ■ 


> 

w 


AMPLITUDE  (dB)  AMPLITUDE  (dB) 

Figure  3.  Amplitude  distribution  histograms  of  events  recorded  at  different  load  levels 
X  100%)  during  quasi-static  loading  to  failure  of  graphite /epoxy 
[0/±45/90]2s  laminates  containing:  a)  non-visual  impact  damage;  and  b)  visual  im¬ 

pact  damage.  The  larger  the  damage  is,  the  larger  the  amount  of  friction  generated 
emission  and  the  higher  the  relative  load  level  at  which  damage  progresses. 


Figure  4.  Event  amplitudes  as  a  function  of  applied  load  recorded  during  quasi-static  loading 
to  failure  of  the  same  specimen  shown  in  Figure  3,  containing:  a)  non-visual  im¬ 

pact  damage;  b)  visual  impact  damage.  In  the  case  of  non-visual  impact  damage  high 
amplitude  events  initiate  at  relatively  lower  load  levels  indicating  early  initia¬ 
tion  of  damage  progression. 
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SPEC.  NO.  HQI-13 
Vi  =  39.3  m/s 
a,  =  273  MPa 
f°  =  1.0  Hz 


1.  (10-100%)a^ 

2.  (10-60%)a^ 

3.  (10-95%) a. 


LOAD  RANGE  =  (10-100%) 
E  =  630  EVENTS 


NUMBER  OF  CYCLES 


LOAD  RANGE  = 

(10-60%) a 
E  =  340  EVENTS 


AMPLITUDE  (dB) 


AMPLITUDE  (dB) 


ENERGY  COUNTS 


ENERGY  COUNTS 


DURATION  [usee] 


DURATION  [ysec] 


COUNTS 

COUNTS 

Figure  5.  Distribution  of  acoustic  emission  source  Intensities  for  events  accumulated  during 
the  initial  part  of  fatigue  loading  (R=0,1)  of  graphlte/epoxy  [  /-  /  12s 

nate:  a)  events  generated  through  the  entire  load  range;  b)  events  generated  at 

the  lower  part  of  the  load  range.  Distinction  between  emission  gene raged  by  frlc 
tion  and  b?  actual  damage  growth  may  require  analysis  of  several  source  Intensities 
and  load  ranges. 


AMPLITUDE  (dB)  DURATION  [psec]  ENERGY  COUNTS  AMPLITUDE  (dB) 


1  Volt  =  45.5  MPa 


SPEC.  NO.  HQI-13  0  =  273  MPa 

V.  =  39.3  m/s  f  =  1.0  Hz 

1 


NUMBER  OF  CYCLES  NUMBER  OF  CYCLES 


NUMBER  OF  CYCLES  NUMBER  OF  CYCLES 


a.  LOAD  RANGE  =  (10-100%)a^ 


b.  LOAD  RANGE  =  (10-60%) a, 

d 


Figure  6.  Acoustic  emission  source  intensities  as  a  function  of  applied  dynamic  stress  and 
number  of  cycles  during  the  initial  part  of  fatigue  loading  of  the  same  specimen 
shown  in  Figure  5;  a)  events  generated  through  the  entire  load  range;  b)  events 
generated  at  the  lower  part  of  the  load  range.  High  source  intensities  occur  only 
at  the  upper  part  of  the  load  range  indicating  damage  accumulation.  At  the  lower 
part  of  the  load  range  all  source  intensities  are  of  medium  and  low  ranges,  all  of 
which  are  generated  by  friction.  These  results  were  obtained  throughout  the  fa¬ 
tigue  loading. 
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RELATIUE 


S.N,  AQI-7  =  34.8  m/s 


S.N,  AQI-15  *  36.0  m/s 


S.N.  AQI-10  =  61.9  m/s 


“7*e  6>0  w 

FREQUENCY  [KHz] 


Figure  7.  Frequency  response  of  impact-damaged  graphite/epoxy 
[0/±45/90]2s  laminate  and  X-radiographs  recorded  for 
three  different  Impact  velocities.  Results  indicate 
that  the  frequency  response  of  the  subject  laminate  is 
sensitive  to  damage  severity. 


inPACT  UELOCITY  Cn^S3 

Figure  8.  Peak  amplitude  of  average  signal  level  (A.S.L.)  of  the  frequency  response  (ex¬ 
amples  of  which  are  shown  in  Figure  7)  measured  at  350  KHz  versus  Impact  velo¬ 
city  (graphlte/epoxy  [0/145/90]2s^  •  peak  amplitude  is  highly  sensitive  to 

impact  damage. 
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Fiber  reinforced  thermoplastic  resin  composites  are  manufactured  by  heating  the  composite  above  the  “melting” 
point  of  the  resin  followed  by  rapid  cooling.  Pressure  is  also  applied  during  cooling.  The  cooling  rate  determines  the 
crystal  structure  of  the  resin.  The  applied  pressure  assures  the  proper  resin  content  and  fiber  distribution.  Hence,  both 
the  cooling  rate  and  the  applied  pressure  affect  strongly  the  m^hanical  properties  of  the  composite.  Therefore,  the 
appropriate  cooling  rates  and  pressures  must  be  employed  to  achieve  the  required  mechanical  properties. 

The  objectives  of  this  investigation  is  to  determine  the  effects  of  the  processing  parameters  on  the  mechanic^  prop¬ 
erties,  and  to  establish  methods  for  selecting  the  proper  processing  variables  (cooling  rate,  pressure)  for  each  application. 
To  achieve  these  objectives,  a  model  is  being  developed  which  simulates  the  manufacturing  process,  and  relates  the 
process  variables  to  the  thermal,  chemical,  and  mechanical  properties.  The  model  consists  of  two  sub-models.  _  The 
“thermo-chemical”  submodel  relates  the  cooling  rate  to  the  crystallinity  and  the  mech^ical  properties  of  the  finished 
product.  The  “fiow”  submodel  relates  the  applied  pressure  and  the  resin  and  fiber  distributions. 

Thus  far,  the  “thermo-chemical”  submodel  has  been  developed.  This  model  provides  relationships  between  the 
cooling  rate  and  crystallinity,  and  between  the  crystallinity  and  mechanical  properties.  The  model  was  implemented 
with  a  “user-friendly”  computer  code  suitable  for  generating  numerical  results  for  thermoplastic  resin  composites  and, 
in  particular,  for  PEEK  resin  composites. 

Tests  were  also  performed  with  PEEK  resin.  In  these  tests,  first  the  relationship  between  the  cooling  rate  and 
crystallinity  was  established  using  differential  scanning  calorimetry.  Second,  the  tensile,  compressive  and  shear  strengths 
and  moduli  were  measured  as  functions  of  the  cooling  rate.  These  data,  together  with  the  computer  code,  can  be  used 
to  determine  the  optimum  cooling  rates  to  be  used  for  PEEK  resin  composites. 
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determine  optimum  process  variables 
(cooling  rate,  pressure) 
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APPROACH 


•  THERMOCHEMICAL  MODEL 
COOLING  RATE 

i 

CRYSTALLINITY 

I 

V 

MECH. PROPERTIES 


O  FLOW  MODEL 


PRESSURE 
RESIN  FLOW 

i 

FIBER-RESIN  DISTR. 

I 

V 

MECH. PROPERTIES 


THERMOCHEMICAL  MODEL 
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HEAT  FLOW 


THERMOCHEMICAL  MODEL 


Absolute  Crystallinity 


MECHANICAL  PROPERTIES 


MECHANICAL  PROPERTIES 


1 .  Tensile  Properties 


a 

ot- 


0  24  0  26  0  28 

Cryslaninily 


2.  Compressive  Properties 


Definiuon  of  Compressive  Sirengih 


♦  Average 
I  Spread 


0  24  0  26  0  :? 

Cryslallinily 


PEEK  150P 


0  24  0  26  0  29 

Crystallinity 


MECHANICAL  PROPERTIES 


3.  Shear  Properties 


SUMMARY 


0  24  0  26  0  28 

Crystallinity 


•  Average 
I  Spread 


PEEK  150P 


0  24  0  26  0  29 

Crystallinity 


PROGRESS  TO  DATE 

1 )  GENERAL  PROCEDURE  FOR  ESTABLISHING 
PROCESS  VARIABLES 

2)  "USER  FRIENDLY  COMPUTER  CODE  FOR 
RELATING  COOLING  RATE, 
CRYSTALLINITY  AND  MECHANICAL 
PROPERTIES 

3)  TEST  METHODS  FOR  DETERMINING 
CRYSTALLINITY 

4)  DATA  AND  EMPIRICAL  CORRELATIONS 
FOR  PEEK  RESIN 

•  CRYSTALLINITY  AS  A  FUNCTION  OF 

COOLING  RATE 

•  MECHANICAL  PROPERTIES  AS  A 

FUNCTION  OF  COOLING  RATE 
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ABSTRACT 

There  have  been  many  studies  of  the  dynamic  behavior  of  composite  materials  and 
structures.  However,  the  majority  of  these  studies  have  been  concerned  with  the 
response  of  laminates  to  localized  impact.  Tolerance  of  the  laminate  to  delamination, 
fiber  breakage,  and  other  forms  of  damage,  and  the  influence  of  any  damage  on  the 
strengths  have  been  the  principal  concerns  of  these  past  studies.  In  contrast,  little 
has  been  done  to  study  large  deformation  transient  dynamic  response  of  structural 
elements  as-a^whole,  that  is  global  response  as  opposed  to  local  response.  However, 
with  the  increasing  usage  of  composite  materials  for  major  structural  elements  in 
aircraft,  more  research  is  needed  to  understand  the  response  of  structural  elements  to 
impacts  which  could  simulate  hard  landings,  wheels-up  landings,  or  other  short  duration 
but  intense  loading  events  that  will  excite  the  entire  structure,  not  just  a  localized 
region.  In  general,  it  is  expected  that  the  response  of  the  structural  elements  to 
these  crash-related  loads  will  involve  large  deformations  and  failure.  In  addition, 
energy  absorption  or  the  lack  of  it  during  failure  is  an  important  variable  if  damage 
to  other  parts  of  the  structure  and  possible  injury  to  occupants  is  to  be  minimized. 
Because  of  these  concerns  with  future  structures,  work  is  in  progress  to  study  the 
large  deformation  response  of  simple  composite  structural  elements  to  intense  dynamic 
loadings. 


A  research  program  has  been  formulated  to  investigate  the  response  characteristics 
of  generic  composite  components  to  simulated  crash  loadings.  This  program  has  been 
arranged  to  focus  on  three  levels:  the  laminate  level  for  material  properties  such  as 
energy  absorbing  qualities  and  the  behavior  of  skin  materials;  the  element  level 
focusing  on  more  complex  geometry  and  behavior  of  beams,  frames  (rings),  arches,  and 
panels;  and  the  substructure  level  dealing  with  cylindrical  shells,  floors,  and  larger- 
scale  components.  Scaling  studies  will  also  be  included. 


The  goal  of  research  on  the  generic  components  is  to  provide  a  data  base  and 
understanding  of  generic  composite  component  behavior  subjected  to  crash  loading 
conditions  supported  by  validated  analytical  methods.  To  help  achieve  this  goal,  in- 
house  research,  contractural  efforts,  and  university  grants  are  included  in  the 
program. 


This  presentation  will  summarize  results  from  three  of  the  research  areas  under 
the  composite  crash  dynamics  program  that  include:  beam  impact  studies,  composite 
fuselage  frame  impact  studies,  and  abrasive  loadings  of  skin  materials. 
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FULL-SCALE  TESTING 

•  Demonstrate  structural  performance  of  composite  full-scale  aircraft  under  crash 
loadings 
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TYPICAL  SKIN  COUPON  ABRASION  TEST  SPECWEN 
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DYNAMICS  AND  AEROELASTICITY  OF  COMPOSITE  STRUCTURES 
John  Dugund ji 


Gun-Shing  Chen 

Technology  Laboratory  for  Advanced  Composites 
Department  of  Aeronautics  and  Astronautics 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139 


ABSTRACT 

In  previous  investigations  at  the  aeroelastic  flutter  and  divergence 

behavior  of  a  series  of  unswept  and  forward  swept,  graphite/epoxy  cantilever  wings 
were  investigated  in  a  small,  low-speed  wind  tunnel.  The  wings  were  six-ply 
graphite/epoxy  plates  and  had  strong  bending^twisting  coupling  (D]_5  terms).  By 
adjusting  the  bending-torsion  coupling,  the  divergence  tendency  of  the  forward  swept, 
cantilever  wings  could  be  eliminated  and  the  flutter  speed  raised  considerably.  See 
Refs.  1  and  2. 

Presently,  an  investigation  is  being  made  into  the  effects  of  rigid  body  aircraft 
modes  on  the  aeroelastic  behavior  of  forward  swept  wings.  It  has  recently  been 
pointed  out  in  Refs.  3  and  4,  that  for  forward  swept  wings,  the  rigid  body  inodes  may 
possibly  couple  with  the  wing  bending  mode  to  cause  a  new  low  frequency  body-freedom- 
flutter."  Accordingly,  a  complete,  two-sided  30°  forward  swept  wing  aircraft  model 
was  constructed  and  mounted  with  low  friction  bearings  in  both  pitching  and  trans¬ 
lation,  inside  the  M.I.T.  low  speed  acoustic  wind  tunnel.  The  wind  tunnel  had  a 
1.5  x  2.3  m  (5  X  7  ft.)  test  section  and  could  reach  velocities  of  30  m/s.  Four 
different  ply  layup  wings  could  be  interchanged  on  the  model,  namely  [02/90)3/ 
[+152/0]c,  [+302/0]s  and  [-152/0]s.  These  wing  surfaces  were  the  same  ones  used  in 

the  previous  cantilever  tests  [Refs.  1  and  2),  and  thus  the  present  tests  complemented 
the  previous  cantilever  tests  and  isolated  the  effects  of  rigid  body  motions. 

The  wind  tunnel  tests  included  measurement  of  the  static  lift  and  moment 
characteristics  (done  at  low  speeds)  and  the  dynamic  stability,  flutter,  and 
divergence  testing  at  higher  speeds.  For  all  free  flying  tests,  the  model  was  set 
to  a  low  trim  angle  of  attack,  and  aircraft  vertical  height.  Also,  TV  movies  were 
taken.  Body-freedom-flutter  was  encountered  for  some  configurations  as  well  as 
torsional  stall  flutter,  bending  stall  flutter,  and  dynamic  instability.  Examples 
are  given  of  various  flutter  and  dynamic  instabilities  encountered.  These  experi¬ 
mental  tests,  along  with  the  previous  cantilever  wing  tests  and  with  corresponding 
analytical  analyses,  should  provide  insight  into  the  actual  aeroelastic  behavior  of 
forward  swept  wing  aircraft  in  free  flight. 
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DYNAMICS  AND  AEROELASTICITT  OF 
COMPOSITE  STRUCTURES 


OBJECTIVES 


John  Dugundji 
Gun-Shing  Chen 


Investigate  the  effects  of  rigid  body  aircraft  f?iodes 
on  the  flutter  and  divergence  of  foroard  swept^ 
graph Ite/epoxy  wings 


Explore  nonlinear  effects  of  large  angle  of  attack  and 
stall  flutter 


•  Obtain  insight  into  actual  aeroelastlc  behavior  of 
forv/ard  swept,  aeroelastically  tailored  aircraft  in  free 
Technology  Laboratory  for  Advanced  Composites  flight 

Department  of  Aeronautics  and  Astronautics 
Massachusetts  Institute  of  Technology 
Cambridge.  Massachusetts  02139 


APPROACH 


Build  complete,  two-sided  30°  forward  swept  wing  aircraft 
model  with  rigid  pitch  and  rigid  translation  copability 


Obtain  experimental  data  on  model  in  low  speed  wind  tunnel, 
and  compare  with  corresponding  cantilever  wing  tests 


•  Perform  analytical  flutter  and  divergence  calculations 
including  effects  of  rigid  body  modes 
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CURRENT  STATUS 


•  Currently^  analyzing  results  of  model  wind  tunnel  tests 

•  Observed  body-freedom-flutter.  torsional  stall  flutter, 
dynamic  instability,  nonlinear  phenomena 

•  Currently,  performing  analytical  flutter  and  divergence 
analyses  to  assess  experimental  results,  compare  with 
cantilever  results,  and  assess  extent  of  nonlinear 
phenomena 


Theoretical  Stability  Plot.  102/9015 
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AN  ELASTIC  STRESS  ANALYSIS  OF  THE  DEBOND 
FRONT  IN  DCB  SPECIMENS 
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K  .  N  .  Shi  vak lima r* 
I .  S.  Raju* 


NASA  Langley  Research  Center 
Hampton,  Virginia  23665-5225 


ABSTRACT 


Double  cantilever  beam  (DCB)  specimens  are  widely  used  to  measure  mode  I  frac¬ 
ture  toughness  of  adhesive  debonding  and  composite  delamination.  Various  combina¬ 
tions  of  adherend  materials  and  specimen  dimensions  are  currently 

these  DCB  specimens.  These  specimen  parameters  can  influence  the  stress  distribu 
tion  ahead  of  the  debond  front  (ref.  1)  and,  therefore,  should  also  ^ 

influence  fracture  toughness.  This  can  be  especially  ' . 
Whereas  brittle  resins  fail  under  the  influence  of  very  localized  stresses  ® 

characterized  by  the  debond  front  st  res  s -i  ntens  i  ty  factor,  tough  resins 
extensively  ahead  of  the  debond  front  before  fracturing.  The  stresses  ahead  of  the 
debond  govern  the  extent  of  yielding  and  the  corresponding  plastic  energy  dissipa¬ 
tion  that  provides  fracture  toughness.  The  objective  of  the  present _ study  was  to 
evaluate  several  parameters  that  may  influence  the  peel  stress  distribution  ahead 
of  the  debond  front. 

Aluminum  and  graphite/epoxy  DCB  specimens  were  studied  using  f i ni te -el ement 
analyses.  The  parameters  that  were  varied  included  adherend  and  adhesive  thick¬ 
nesses  as  well  as  adherend  longitudinal  and  transverse  (thi ckness -di recti  on )  stiff¬ 
nesses.  Results  were  plotted  as  a  (peel)  stress  distributions  ahead  of  the 

y 

debond  front . 

Neither  the  adherend  thickness  nor  the  longitudinal  stiffness  had  f^uch  influ- 
ence  on  the  peel  stress  distributions,  despite  their  strong  influence  on  the  stress 
intensity  factors.  In  contrast,  both  adherend  transverse  stiffness  and  adhesive 
thickness  had  a  noticeable  influence  on  the  stress  di s t n bu t i 9ns ,  but  did  not 
affect  stress  intensity  factors.  Plastic  zone  sizes  were  estimated  for  a  wide 
range  of  adhesive  thicknesses  and  were  found  to  have  a  peak  ^al ue  wi thi n  thi s 
range.  This  trend  agrees  qualitatively  with  the  well  known  dependence  that  DCB 
fracture  toughness  has  on  adhesive  thickness. 
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STRESS  DISTRIBUTION  AHEAD  OF  DELAMINATION 


COMPARISON  OF  STRESS  DISTRIBUTIONS 


ANALYSIS  OF  STRESS  DISTRIBUTION  FOR  ALUMINUM  DCB 


Beam  on  elastic  foundation 

Finite  element  onolysls 


lU 

10“^^  10“^  lO"**  10”^  10“^  10“'  lO'" 

Xy  mu 


d)  t  =  0.36  m 


e)  t  =  0.86  Compact  tension 


YIELD  ZONE  SIZES 


Sim"ARY 


0  ADHESIVE  Cy  (PEEL)  STRESSES  HERE  HIGHER  THAN 

THOSE  FOR  AONOLITHIC  CASE. 

0  ADHEREND  FLEXURAL  STIFFNESS  HAD  LITTLE  EFFECT 
ON  (Ty  DISTRIBUTION. 

0  ADHERENO  TRANSVERSE  STIFFNESS  (IN  THICKNESS 
DIRECTION)  HAD  SIGNIFICANT  EFFECT  ON  CT^ 

DISTRIBUTION 

0  PLASTIC  ZONE  SIZE  VARIED  WITH  ADHESIVE  THICKNESS 
AND  HAD  A  PEAK  VALUE. 


Adhesive  thickness^  t,  mn 


STATISTICAL  LIFE  PREDICTION  OF  COMPOSITE  STRUCTURES 
D.  Jones  and  J.  Yang 
The  George  Washington  University 

Material  not  received  in  time  for  publication. 
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ABSTRACT 


Delamination  is  the  most  commonly  observed  damage  mode  in  composite  lami¬ 
nates.  In  thin  plates  under  low-velocity  impact  from  point  loads,  matrix  cracks 
will  form  first  in  the  laminate  plies,  and  delaminations  will  then  grow  from  these 
matrix  cracks  where  they  terminate  at  ply  interfaces.  In  this  study,  delamination 
growth  from  matrix  cracks  in  laminates  subjected  to  bending  loads  was  evaluated. 
Using  laminate  plate  theory,  simple  eguations  were  derived  which  relatethe  strain 
energy  release  rate  associated  with  the  delamination  growth  to  the  applied  load  and 
laminate  stiffness  properties.  Several  bending  problems  and  boundary  conditions 
were  considered.  In  order  to  isolate  the  matrix  crack  and  model  the  problem  2- 
dimensionally,  different  layups  of  graphite-epoxy  (5208/T300)  laminates  were 
inodsled  and  tested.  All  tested  layups  were  symmetric  and  had  90  plies  on  the 
outside.  For  such  a  laminate  under  a  bending  load,  matrix  cracks  will  form 
parallel  to  the  fibers,  through  the  90®  plies,  on  the  tension  surface  of  the  lami¬ 
nate.  Delaminations  then  form  and  grow  in  the  interface  between  the  Inmost  cracked 
90°  ply  and  the  adjacent  ply.  To  simulate  pure  bending,  specimens  were  loaded 
across  their  width  by  pins  which  were  free  to  rotate  on  ball  bearings,  in  a  4-point 
bending  apparatus.  A  three-point  bend  test  was  performed  in  the  same  type  appa¬ 
ratus  by  supporting  the  specimen  on  two  pins  and  loading  it  at  the  midspan.  To 
simulate  the  deformations  in  laminates  subjected  to  low-velocity  impact,  specimens 
were  clamped  at  the  ends  and  loaded  across  their  width  at  the  midspan.  For  all 
configurations,  the  laminates  were  loaded  until  delamination  was  first  observed. 

The  delamination  loads  predicted  by  the  analysis  correlated  well  with  the  measured 
delamination  loads  for  the  different  layups  and  boundary  conditions, 

*PRC  Kentron,  I'nc. ,  Hampton,  Virgi  nTa" 
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ABSTRACT 


The  mechanics  of  transverse  cracking  in  an  elastic  fibrous  composite  ply  is  explored,  first  for 
low  crack  density,  and  then  during  progressive  cracking  under  increasing  load.  Cracks  are  assumed  to 
initiate  from  a  nucleus  created  by  coalescence  of  localiz^  fiber  debonding 

Conditions  for  onset  of  unstable  cracking  are  evaluated  with  regard  to  it  L  found 

adjacent  plies  of  different  elastic  properties.  Interfacial  damage  is  considered  as  well,  I 
that  crack  propagation  may  take  place  in  two  directions  on  planes  which  are  parallel  to  the  ^ 

perpendicular  to  the  midplane  of  the  ply,  and  that  the  actual  propagation  direction  and  ply  strength 
depend  on  ply  thickness. 

Thin  dies  are  defined  by  the  requirement  that  the  crack  nucleus  extends  across  the 
thicknels  before  Snseiof  unstable  crack  propagation.  Of  course,  plies  with  preexisting  through- the 
thickness  cracks  caused  by  other  damage  modes,  such  as  impact,  should  also  be  ^ ’f 

regardless  of  their  actual  thickness.  Strength  of  thin  plies  is  found  to  be  ^"trol  ed  by  onset  of 
crack  propagation  in  the  fiber  direction;  these  longitudinal  cracks  are  called  type  L  cracks.  The 
strain^energy  release  rate  G(L)  for  these  cracks  is  derived.  This  leads  to  expressions  which  relate 
ply  strength  to  ply  thickness  and  toughness  Gc(L).  Both  simple  and  mixed  loading  modes 
Only  two  strain  states,  normal  tension  and  longitudinal  shear  affect  p  y  failure.  Traisverse  hear 
seems  to  be  neutralized  by  interlocking  asperities  created  on  the  crack  surface  as  the  crack  finds  its 
way  between  closely  spaced  fibers. 

Thick  olies  are  defined  by  the  condition  that  the  width  of  a  crack  nucleus  of  critical  size  is 
much  smaller  than  ply  thickness.  Thus  there  is  no  interaction  between  the  nucleus  and  adjacent  plies. 
The  energy  release^rates  for  the  nucleus  are  derived.  In  addition  to  G(L)  for  type  L  cracks,  we  also 
evaluate  G(T)  that  corresponds  to  type  T  cracks  which  propagate  in  the  transverse  direction, 
perpendicular  to  the  fiber  axis,  across  the  thickness  of  the  ply.  It  is  J' 

there  is  a  strong  preference  for  the  crack  nucleus  to  propagate  first  as  type  ^  P  ^ 

thickness,  and  tLn  change  direction  and  continue  as  type  L  crack  across  the  entire  oaded  area  of  the 
ply.  It  is  recognized  that  the  two  different  crack  directions  involve  different  values  of  ply 
toughness,  G^,  and  that  Gc(L)  <  Gc(T).  However,  experimental  evidence  appears  to  indicate  that 
Gc(T)  <  2Gr(L).  hence  the  strength  of  thick  plies  is  controlled  by  type  T  cracking,  ^ 

strength  o?  thin  plies  is  controlled  by  type  L  cracking.  Expressions  for  thick  ply  . 

obtained  in  terms  of  critical  size  of  crack  nucleus  and  ply  toughness  G(.(T).  Since  the  nucleus  size  is 
not  usually  known,  one  cannot  predict  the  strength  of  thick  plies.  However,  it  is  possible  to  show 
that  this  strength  does  not  depend  on  ply  thickness  for  a  given  history  of  loading  to  failure. 


Prepared  for  the  Eleventh  Annual  Mechanics  of  Composites  Review,  Dayton,  October  22-24.  1985 
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These  results  for  thin  and  thick  plies  compare  well  with  experimentally  measured  dependence  of  ply 
strength  on  ply  thickness.  One  of  the  conclusions  that  can  be  drawn  from  these  results  is  that 
strength  of  all  plies  should  be  evaluated  from  thin  ply  formulae,  to  guard  against  strength  reduction 
In  thick  plies  which  may  be  caused  by  preexisting  cracks. 

In  the  second  part  of  the  presentation,  we  analyse  progressive  cracking  in  a  ply  under  Increasing 
applied  strain.  This  process  Is  essentially  a  continuous  repetition  of  first  ply  failure,  but  both 
formation  of  the  Initial  crack  nuclei,  and  energy  release  rates  G(L)  and  6(T)  for  unstable  cracks  are 
influenced  by  Interaction  between  the  new  crack  and  the  cracks  which  are  already  present.  This  type  of 
crack  Interaction  Is  analyzed  by  several  methods;  the  self-consistent  method,  as  well  as  elasticity 
solutions  for  a  row  of  parallel  cracks  are  employed.  It  Is  found  that  very  similar  crack  Interaction 
effects  are  predicted  by  these  different  approaches.  Of  particular  interest  is  the  evaluation  of  G(L) 
and  G(T)  of  a  crack  which  Is  about  to  form  In  a  ply  that  already  contains  a  certain  density  of  similar 
cracks.  Also,  the  strain  in  the  ligaments  between  existing  cracks  Is  evaluated. 
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STRENGTH  OF  THICK  PLIES  IS  CONTROLLED  BY  TRANSVERSE 
(type  T)  CRACKS  WHICH  ARE  AUTOMATICALLY  FOLLOWED  BY 
LONGITUDINAL  (type  L)  CRACKS. 


PART  1  -  MECHANICS  OF  FIRST  PLY  FAILURE  CRACK  NUCLEATION 
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STRENGTH  OF  THIN  PLIES  -  1  STRENGTH  OF  THIN  PLIES 
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In  the  presence  of  interface  damage: 


STRENGTH  OF  THICK  PLIES 
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INTERACTION  ENERGY  FOR  ROW  OF  PARALLEL  CRACKS 


57 


at  e  =  0:  G°(L)  =  W^/Za 


009 


c 

o 


u 


c 


03 

to 

(tJ 


<u 

s- 


58 


INITIATION  OF  FREE-EDGE  DELAMINATION  IN 
COMPOSITE  LAMINATE  -  PREDICTION  AND  EXPERIMENT 

Ran  Y.  Kim 

University  of  Dayton 
Research  Institute 
Dayton,  Ohio  45469 

and 
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ABSTRACT 

Delamination  along  the  straight  free  edge  of  composite  laminates  under  an  inplane 
uniaxial  load  has  been  observed  since  the  early  1970 *s.  The  growth  of  this  delamina¬ 
tion  under  load  often  results  in  a  reduction  of  strength  and  stiffness  of  the  lami¬ 
nate.  Since  then  a  large  number  of  papers  have  reported  on  the  free-edge  problem  in 
composite  laminates,  indicating  that  free-edge  delamination  is  attributed  to  the 
existence  of  interlaminar  stresses  which  are  highly  localized  in  the  neighborhood  of 
free  edge  under  an  inplane  loading.  As  a  result  of  this  analytical  work  [1,2],  the 
nature  of  interlaminar  stresses  with  regard  to  magnitude  and  sign  of  each  interlami¬ 
nar  stress  component  is  reasonably  well  understood.  In  spite  of  the  existence  of  the 
analytical  models  that  provide  interlaminar  stress  distributions  in  the  free-edge 
region,  no  successful  prediction  method  for  the  onset  of  delamination  was  reported. 

In  the  present  work  an  effort  has  been  made  to  adequately  predict  the  onset  of 
delamination  using  stress  analysis  in  conjunction  with  a  failure  theory. 

The  first  phase  of  this  work  was  to  determine  the  interlaminar  normal  stress  (a  ) 
experimentally.  Strain  in  the  thickness  direction  was  measured  at  the  free  edge  by 
employing  a  miniature  strain  gage,  and  Cz  is  calculated  using  three-dimensional  con¬ 
stitutive  relations  in  conjunction  with  axial  and  transverse  strain  [3].  The  experi¬ 
mental  results  compare  very  well  with  the  analytical  results  calculated  from  the 
global-local  model  developed  by  Pagano  and  Soni  [4] . 

The  second  phase  of  this  work  was  to  develop  a  methodology  to  adequately  predict 
the  onset  of  delamination  due  to  either  interlaminar  normal  stress  [5]  or  shear  stress 
alone  [6] .  The  global-local  laminate  variational  model  has  been  used  to  calculate  the 
interlaminar  stresses,  and  the  average  stress  components  over  one  layer  from  the  free 
edge  are  assumed  to  be  the  effective  stress  level  acting  at  the  free  edge.  The  pre¬ 
dicted  values  in  most  cases  are  fairly  well  compared  with  the  experimental  results  for 
a  variety  of  laminates,  and  the  predicted  interface  of  delamination  is  the  same  as  the 
experimentally  observed  one.  The  quadratic  failure  theory  [7]  has  been  utilized  to 
improve  the  prediction  capability  in  the  case  where  both  normal  and  shear  stress 
components  are  significant  [8] . 
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OBJECTIVE 

TO  DEVELOP  A  METHODOLOGY  FOR  PREDICTING  THE  FREE-EDGE  DELAMINATION  THRESHOLD 
STRESS  AND  LOCATION  IN  COMPOSITE  LAMINATES  UNDER  APPLIED  UNIAXIAL  LOAD 


APPROACH 

CALCULATE  INTERLAMINAR  STRESS  COMPONENTS  AT  FREE  EDGE  BY  GLOBAL-LOCAL  MODEL 
OBSERVE  DELAMINATION  EXPERIMENTALLY 
DETERMINE  <7^ 

DETECT  ONSET  OF  DELAMINATION 
EXAMINE  THE  FRACTURED  SURFACE 

EMPLOY  AVERAGE  STRESS  FAILURE  THEORY 
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STRESS  AND  STRENGTH  ANALYSIS  AT  DELAMINATION 
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QUADRATIC  POLYNOMIAL  FAILURE  CRITERION  STRENGTH  RATIO  R  AT  ONSET  OF  DELAMINATION 

^  USING  QUADRATIC  FAILURE  CRITERION 
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DELAMINATION  DUE  TOOL  UNDER 
APPLIED  COMPRESSION 


DELAMINATION  DUE  TO  T  UNDER 
APPLIED  COMPRESSION 


CONCLUSIONS 


THE  EXPERIMENTALLY  DETERMINED  Oz  CORRELATES  CLOSELY  WITH  THE  ANALYTICAL 
RESULT. 

THE  PREDICTED  ONSET  OF  DELAMINATION  LOADS  EMPLOYING  AVERAGE  STRESS 
FAILURE  CRITERIA  ARE  CLOSE  TO  THE  OBSERVED  VALUES  WHERE  SINGLE 
PREDOMINANT  INTERLAMINAR  STRESS  COMPONENT  EXISTS. 

QUADRATIC  FAILURE  CRITERION  PREDICTS  STRENGTH  RATIOS  CLOSE  TO 
OBSERVED  ONES  FOR  A  NUMBER  OF  LAMINATES. 

THE  MODE  AND  LOCATION  OF  DELAMINATION  IN  COMPOSITE  LAMINATES  ARE 
PREDiaED  ACCURATELY.  THE  PREDICTED  MODES  OF  DELAMINATION  DUE  TO 
OR  Tx2  can  be  clearly  IDENTIFIED  BY  SEM. 


FURTHER  WORK  IS  REQUIRED  TO  INVESTIGATE  THE  EFFECT  OF  ORIENTATION  ON 
INTERLAMINAR  SHEAR  STRENGTH. 


STRENGTH,  DEFLECTIONS  AND  IMPACT  DAMAGE  IN  ADVANCED  COMPOSITE  STRUCTURES 

Paul  A.  Lagace 

TECHNOLOGY  LABORATORY  FOR  ADVANCED  COMPOSITES 
Department  of  Aeronautics  and  Astronautics 
Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139 


ABSTRACT 

The  plans  for  and  the  progress  on  two  three-year  contractual  efforts 
concerning  composite  materials  will  be  presented.  Although  the  two  major  work 
areas  are  related,  they  will  be  treated,  for  the  purpose  of  presentation, 
separately.  The  first  area  involves  the  "Strength  and  Deflections  of  Advanced 
Composite  Sandwich  Structures".  The  second  area  is  entitled  "The  Sensitivity 
of  Kevlar/Epoxy  and  Graphite/Epoxy  Structures  to  Damage  from  Fragment  Impact". 


Sandwich  construction  has  been  used  in  the  aerospace  industry  for  a  number 
of  years.  The  current  effort  seeks  to  investigate  such  construction  for 
application  to  high  efficiency,  ultralight  aircraft.  This  involves  the 
investigation  of  sandwich  structures  with  relatively  thin  graphite/epoxy  face 
sheets.  Due  to  this  construction,  a  number  of  failure  modes  must  be 
considered:  ply  or  sublaminate  buckling,  skin  debonding,  core  failure,  plate 
buckling,  and  facesheet  fracture.  This  three-year  program  has  thrusts  in  three 
principal  areas;  the  development  of  analytical  tools  necessary  to  predict  the 
strength  and  deflection  of  graphite/epoxy  sandwich  structures;  experimentation 
on  sandwich  plates  to  determine  their  buckling  and  deflection  characteristics 
to  compare  with  the  analysis;  and  experimentation  on  sandwich  plates  to 
determine  their  strength  properties. 

The  work  during  the  last  nine  months  has  centered  on  the  development  of 
the  analytical  methodology.  The  analysis  is  based  on  a  Rayleigh-Ritz  type 
discretization  of  the  displacement  fields  and  uses  direct  minimization  of  the 
potential  energy.  All  eighteen  terms  in  the  A,  B,  and  D  matrices  are  taken 
into  account  so  that  anisotropic  plates  can  be  analyzed.  Mindlin  plate  theory 
is  used  to  incorporate  the  effects  of  transverse  shear  and  the  effect  of 
initial  imperfections  can  also  be  evaluated  as  the  program  uses  Marguerre's 
shallow  shells  theory.  The  loading  on  the  plate  can  be  longitudinal 
compression,  transverse  compression,  shear,  or  any  combination  thereof.  The 
program  solves  the  general  buckling/postbuckling  problem  via  an  optimal  search 
and  a  typical  plate  takes  about  20  seconds  of  CPU  time  on  a  VAX  11/782  computer 
for  analysis.  Typical  results  are  presented  and,  where  available,  agree  well 
with  other  reported  analytical  results. 

An  experimental  program  has  also  been  devised.  The  purpose  of  the  program 
is  to  investigate  the  effects  of  various  core  materials  on  the  deflections  and, 
ultimately,  strength  behavior.  The  experimental  technique  is  discussed  and 
preliminary  results  presented  which  correlate  well  with  the  analysis  for  these 
configurations . 


The  response  of  composite  laminates  to  impact  is  an  important 
consideration  in  assessing  the  damage  tolerance  of  a  composite  structure. 
Composites  are  sensitive  to  impact  due  to  their  tendency  to  delaminate . 
Furthermore,  this  damage  often  goes  undetected  although  it  may  cause 
considerable  strength  reduction  in  the  composite  part.  The  effort  in  this 
three-year  program  is  directed  toward  first  establishing  the  basic  response  of 
composites  to  impact  via  generic  specimens  (i.e.  coupon  type),  analysis  of  the 
impact  event,  and  the  response  of  the  composite  to  the  damage  induced  by 
impact.  Once  the  basic  mechanisms  are  established  and  better  understood,  the 
work  will  progress  to  structures  typically  used  in  aircraft  such  as  stiffened 
panels,  pressurized  cylinders  (which  model  fuselages)  and  the  like. 
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The  analysis  of  the  impact  event  has  been  broken  down  into  parts.  The 
treatment  involves  either  analysis  or  experimentally  determined  data  for  the 
contaT  law  and  a  plate  solution  for  the  overall  parameters  in  order  to  model 
the  local  contact  problem  of  a  static  indentor.  This  is  used  with  a  plate 
solution  to  determine  the  dynamic  response  on  the  plate  to  the  h 
These  two  analyses  yield  the  stress  and  strain  fields  which  are  then  used  in 
coniunction  with  failure  ■  criteria  to  determine  the  damage  caused  the 
impact.  A  modified  Bessel  function  approach  is  used  to  solve  the  axisymmetric 
contact  problem.  The  analytical  solution  shows  excellent  cor relation^wi th 
exoe^iLntal  data  presented  in  the  literature.  The  dynamic  problem  is  solved 
via  a  Rayleigh-Ritz  formulation  and  includes  shear  deformation.  The  technique 
is  currently  being  implemented. 


An  extensive  experimental  program  is  underway  as  well.  An  impact  gun  has 
been  Ltuf  along  with  monitoring  equipment  and  this  will  be  discussed  A 
numCer  of  graphite/epoxy  coupons  have  been  impacted  and  then  sectioned  to 
characterize  the  damage.  Typical  damage  modes  will  be  shown.  Some  tests  to 
failure  have  also  been  conducted  on  impacted  specimens  and  these  preliminary 

results  will  also  be  discussed.  Additionally,  work  has  been  accomplished  on 

modelling  damage  due  to  impact  via  delaminations  implanted  via  ^o£ 

r  number  of  other  authors  have  looked  at  the  compressive  response  of 
qraphite/epoxy  with  implanted  delaminations.  in  this  investigation,  t 

initial  wo?k  has  been  conducted  under  tension  in  order  to  separate  out  the 

basic  effect  the  delamination  has  on  the  fracture  ^ 

structural  effect  the  delamination  induces  due  to  local  ply/sublaminate 
buckling.  The  preliminary  results,  which  will  be  presented,  indicate  that  a 
SJaiiSa^ion  located  at  th^  midplane  of  a  is  benign  in  that  no 

strength  and  little  stiffness  changes  are  noted.  However,  the  cases  wne 

delaminations  are  implanted  at  every  ply  interface,  which  is 
the  damage  in  an  impacted  specimen,  show  strength  reductions  of  30%  for  a  30 
circular  delamination  in  a  70  mm  wide  specimen. 


in 

research 


addition  to  the  results  attained 
will  be  outlined  in  each  case. 


to  date, 


the  next  steps  in  the 
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ANALYSIS  OF  COMPOSITE  LAMINATES 
WITH  DELAMINATIONS  UNDER  COMPRESSION 
LOADING 


HAN-PIN  KAN 


NORTHROP  CORPORATION 
AIRCRAFT  DIVISION 
HAWTHORNE,  CALIFORNIA  90250 


ABSTRACT 


Delaminations  in  composite  structures  are  of  great  concern  to  aircraft  designers  and  analysts. 

This  is  because  of  the  possible  local  buckling  and  growth  of  the  delaminations  when  such  a  structure 
is  subjected  to  in-service  loading.  An  extensive  technology  assessment  conducted  under  Reference  1  has 
identified  that  impact  damage  and  delamination  are  the  most  critical  damage  modes  that  degrade  struc¬ 
tural  strength  and  fatigue  life.  Unlike  metallic  structures,  where  the  material  defects  propagate 
under  tensile  loading,  delaminations  in  composite  laminates  usually  do  not  propagate  under  tension. 
However,  under  compressive  static  and  fatigue  loading,  local  buckling  of  the  delaminated  region  may 
occur  below  the  compression  strength  of  the  laminate.  The  local  buckling  of  the  delaminated  region 
provides  a  means  of  releasing  strain  energy  and  thus  making  energy  available  for  the  delamination  to 
grow.  To  design  composite  structures  for  damage  tolerance,  analytical  methods  are  needed  to  predict 
the  initial  buckling,  to  compute  the  strain-energy-release  rate  and  to  determine  the  mode  of  failure 
of  the  structure. 

In  the  present  paper,  an  analytical  method  to  determine  the  initial  buckling  of  the  delaminated 
region  and  to  compute  the  strain-energy-release  rate  is  presented  and  the  interaction  of  laminate 
failure  mode  is  discussed.  Three  types  of  delaminations,  namely,  through-the-^idth,  circular  and 
elliptical  delaminations  are  considered.  The  influence  of  delamination  size  and  depth  on  the  initial 
buckling  and  strain-energy -release  rate  is  investigated. 

The  initial  buckling  analysis  is  based  on  the  energy  method  employing  the  Rayleigh-Ritz  method. 

The  delaminated  laminate  is  modelled  as  a  homogeneous,  elastic,  orthotropic  plate.  The  plate  is 
separated  into  two  layers  in  the  delaminated  region.  The  strain  energy  expression  is  taken  from 
Reference  2,  Out-of-plane  displacement  functions,  satisfying  the  boundary  conditions,  are  assumed 
for  different  delamination  types.  The  initial  buckling  load  is  determined  by  performing  the  necessary 
integrations  and  variations  to  minimize  the  total  potential  energy  of  the  system.  Results  of  the 
analysis  show  that  the  initial  buckling  load  depends  on  the  delamination  size,  shape,  depth  and  the 
laminate  stacking  reference. 

As  the  applied  compression  load  exceeds  the  initial  buckling  load  of  a  laminate  with  a  delamina¬ 
tion,  strain  energy  will  be  released  for  delamination  growth.  The  total  strain-energy-release  rate, 

G,  is  computed.  The  value  of  G  is  obtained  under  a  constant  load  and  according  to  the  usual  defini¬ 
tion  that  G  is  the  energy  required  to  create  a  unit  new  free  surface  in  the  system. 

The  two  layers  of  a  composite  laminate  separated  by  a  delamination  behave  as  a  single  undamaged 
laminate  when  the  applied  load  is  below  the  initial  buckling  load;  hence  the  strain  energy  release 
rate  G  =  0  in  the  prebuckling  range.  Beyond  initial  buckling,  the  value  of  G  depends  on  the  initial 
buckling  load  which  is  a  function  of  the  delamination  size,  shape  and  depth.  The  results  indicate 
that  the  value  of  G,  under  a  constant  compression  loading,  is  very  sensitive  to  delamination  size, 
when  the  delamination  is  slightly  larger  than  the  critical  size  corresponding  to  the  applied  load. 

As  the  delamination  becomes  large  as  compared  to  the  critical  size,  the  value  of  G  approaches  a 
constant. 

The  initial  buckling  of  the  delaminated  region  in  a  laminate  with  delamination  usually  does  not 
cause  failure  of  the  laminate.  The  static  compression  strength  of  a  composite  structure  with  delamina¬ 
tions  depends  on  the  initial  buckling  strength,  the  material  compression  strength  and  interlaminar 
toughness,  and  the  structural  geometries.  The  mode  of  failure  may  be  gross  compression,  local  or 
global  buckling  or  catastrophic  delamination  growth.  A  failure  prediction  technique  is  developed 
which  takes  into  consideration  of  the  failure  mode  interactions.  In  this  technique,  the  failure 
strength  for  different  failure  mode  is  plotted  as  a  function  of  the  delamination  size.  The  laminate 
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failure  strength  is  then  given  by  the  minimum  strength  among  all  the  failure  modes.  By  doing  this 
the  failure  is  divided  into  three  or  more  regions  depending  on  the  delamination  size. 
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ANALYSIS  OF  COMPOSITE  LAMINATES  WITH 


DELAMINATIONS  UNDER  COMPRESSION  LOADING 


HAN-PIN  KAN 


NORTHROP  CORPORATION 
AIRCRAFT  DIVISION 
HAWTHORNE.  CALIFORNIA  90250 


OBJECTIVE 


TO  PREDICT  RESIDUAL  STRENGTH  AND  FAILURE  MODE 
OF  COMPOSITE  LAMINATES  WITH  DELAMINATIONS 
UNDER  COMPRESSION  LOADING 


APPROACH 


•  DETERMINE  INITIAL  BUCKLING  OF  DELAMINATED  REGION 

•  COMPUTE  STRAIN-ENERGY-RELEASE  RATE  FOR 
DELAMINATION  GROWTH 

•  CONSIDER  FAILURE  MODE  FOR  STRENGTH  PREDICTION 
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ANALYSIS  FEATURES 


•  SINGLE  DELAMINATION 


•  ANY  DELAMINATION  SIZE,  SHAPE,  AND  DEPTH 


•  BOUNDARY  VARIATION  OF  STRAIN-ENERGY-RELEASE 
RATE  INCLUDED 

•  MULTIPLE  FAILURE  MODES  FOR  STRENGTH  PREDICTION 


BUCKLING  ANALYSIS 


•  DELAMINATED  LAMINATE  MODELED  AS  TWO 
SEPARATED  LAYERS  OF  HOMOGENEOUS, 
ELASTIC,  ORTHOTROPIC  PLATES 

•  ENERGY  FORMULATION  EMPLOYED 

•  INITIAL  BUCKLING  STRENGTH  OF  THE 
WEAKER  LAYER  IS  OBTAINED  USING 
RAYLEIGH-RITZ  METHOD 

•  INITIAL  BUCKLING  LOAD  DEPENDS  ON  THE 
DELAMINATION  SHAPE,  SIZE,  LOCATION,  AND 
THE  LAMINATE  STACkiNG  SEQUENCE 


THROUGH-THE-WIDTH  DELAMINATION 


LAMINATE  WITH  AN  ELLIPTICAL 
DELAMINATION 


Y 


DISPLACEMENT  FUNCTIONS 


THROUGH-THE-WIDTH  DELAMINATION 
w  =  w(x)  =  Ae®’*  +  +  Csin  a  x  +  Dcos  a  x 

ELLIPTICAL  DELAMINATION 

w  =  w(x,y)  =  — - —  (a*b2  -  b*  x*  -  a*  y*)* 


INITIAL  BUCKLING  LOADS 


THROUGH-THE-WIDTH  DELAMINATION 

Ncr,  =  k(-^)*  •  Di,  .  k=  1.0306 

for  lowest  buckling  mode 

ELLIPTICAL  DELAMINATION 

Ncri  =  +  ^(-|)*  D„] 

Djj  ARE  PLATE  RIGIDITIES  OF  THE  DELAMINATION  LAYER 
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DISPLACEMENT,  (inch 


INITIAL  BUCKLING  STRENGTH 


DELAMINATION  LENSTH-IN. 


DISPLACEMENT  PROFILE 
OF  THE  DELAMINATED  LAYER 


-0.4 


-0.3  -0.2  -0.1  0.0  0.1  0.2  0.3 

DISTANCE  FROM  CENTERLINE  OF  DELAMINATION,  X  (Inch) 

STRAIN-ENERGY  RELEASE 
RATE  ANALYSIS 


0.4 


•  ASSUME  NO  ENERGY  RELEASED  AT  APPLIED  LOAD 
LEVEL  BELOW  INITIAL  BUCKLING  STRENGTH  OF  THE 
DELAMINATED  LAYER 

•  OVERALL  STRAIN-ENERGY- RELEASED-RATE  IS 
DETERMINED  UNDER  CONSTANT  LOAD  CONDITION 

•  STRAIN-ENERGY- RELEASE- RATE 

^  _  acw-u) 

•  G  IS  A  FUNCTION  OF  DELAMINATION  SHAPE,  SIZE, 
LOCATION,  AND  THE  INITIAL  BUCKLING 

LOAD 

•  FOR  CIRCULAR  AND  ELLIPTICAL  DELAMINATIONS, 

G  ALSO  VARIES  AROUND  THE  DELAMINATION 
FRONT 
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STRAIN-ENERGY  RELEASE  RATE 

1(i-PLY  LAMtNATE  3-PLY-OEEP  DELAMINATION 


EFFECTS  OF  STACKING  SECH^NCE  ON 
STRAIN  ENERGY  RELEASE  RATE 


HALF  DELAMINATION  LENGTH  -  IN. 


77 


EFFECTS  OF  PLY  THICKNESS  ON 
STRAtN-ENERGY-RELEASE  RATE 


STRAIN-ENERGY-RELEASE  RATE 
CIRCULAR  DELAMINATION 


78 


KRIPHERAL  VAmTION  OF  STRAIN-ENERGY- 
ROEASE-RATE  CIRCULAR  KLAMIHATION 


EFFECTS  OF  LAMINATE  THICKNESS 
ON  STRAIN-ENERGY-RELEASE-RATE 


EFFECTS  OF  DELAMINATION  DEPTH 
ON  STRAIN-ENERGY-RELEASE-RATE 
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STRAIN-ENERGr-RELEASE-RATE 
ELLIPTICAL  DELAMINATION 


EFFECT  OF  ELLIPTICAL  ASPECT  RATIO 
ON  THE  STRAIN-ENERCT-RELEASE-RATE 


STATIC  FAILURE  MODES  FOR 
LAMINATE  WITH  DELAMINATIONS 


MODE 

DESCRIPTION 

1 

GROSS  COMPRESSION  FAILURE 

2 

INITIAL  BUCKLING  OF  THE  DELAMINATION  FOLLOWED  BY 
COMPRESSION  FAILURE  OF  THE  REMAINING  LAYER 

3 

GLOBAL  BUCKLING 

4 

INITIAL  BUCKLING  OF  THE  DELAMINATION  FOLLOWED  BY 
GLOBAL  BUCKLING  OF  THE  LAMINATE  GLOBAL 

5 

LOCAL  BUCKLING 

INITIAL  BUCKLING  OF  THE  DELAMINATION  FOLLOWED  BY 
LOCAL  BUCKLING  OF  THE  REMAINING  LAYER 

6 

DELAMINATION  GROWTH 

INITIAL  BUCKLING  OF  THE  DELAMINATION  FOLLOWED  BY 
DELAMINATION  GROWTH  AND  EVENTUALLY  FAILED  BY 

MODE  2.  4  OR  5 
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GROSS  COMPRESSION  FAILURE 
(FAILURE  MODE  1) 


INITIAL  BUCKLING  OF  THE  DELAMINATICIN 
FOLLOWED  BY  COMPRESSION  FAILURE 
(FAILURE  MODE  2) 


GLOBAL  BUCKLING 
(FAILURE  MODE  3) 


INITIAL  DELAMINATION 


INITIAL  BUCKLING  OF  THE  DELAMINATION 
FOLLOWED  BY  GLOBAL  BUCKLING 
(FAILURE  MODE  4) 


LOCAL  BUCKLING 
(FAILURE  MOPE  5) 


INITIAL  BUCKLING  OF  THE  DELAMINATION 
FOLLOWED  BY  DELAMINATION  GROWTH 
(FAILURE  MODE  6) 
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FAILURE  MODE  INTERACTION 


SUMMARY 


•  AN  ANALYTICAL  PROCEDURE  HAS  BEEN  DEVELOP, ED  FOR  STATIC 
COMPRESSION  STRENGTH  PREDICTION  OF  COMPOSITE  LAMINATES 
WITH  DELAMINATIONS 

•  the  influence  of  DELAMINATION  CONFIGURATION  ON  THE  FAILURE 
STRENGTH  AND  FAILURE  MODE  HAS  BEEN  INVESTIGATED 

•  THE  INFLUENCE  OF  LAMINATE  THICKNESS.  PLY  THICKNESS  AND 
STACKING  SEQUENCE  ON  STRAIN-ENERGY-RELEASE  RATE  HAS  BEEN 
EXAMINED 
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SIMPLE  RECTANGULAR  ELEMENT  FOR 
ANALYSIS  OF  LAMINATED  COMPOSITE  MATERIALS 

John  D.  Whitcomb 

National  Aeronautics  and  Space  Administration 
Langley  Research  Center 
Hampton,  Virginia  23665-5225 


ABSTRACT 

The  development  of  an  appropriate  finite  element  mesh  is  a  key  step  in  sue- 
cessful  finite  element  analysis.  For  homogeneous  materials  the  mesh  refinement  is 
dictated  by  geometrical  considerations.  The  shape  of  a  structure  should  be  faith- 
fully  modeled.  Also,  extra  mesh  refinement  is  required  in  regions  with  strong 
strain  gradients  caused  by  holes,  cracks,  or  boundary  conditions.  For  laminated 
materials  the  analyst  must  also  account  for  the  different  material  properties  of 
the  various  laminae.  These  ideas  are  illustrated  by  the  laminated  composite  beam 
shown  in  the  first  figure.  Geometrical  considerations  require  very  few  elements 
except  close  to  the  points  where  load  is  applied,  where  strain  gradients  are 
large.  But  since  standard  finite  elements  cannot  account  for  stacking  sequence 
effects,  such  elements  should  not  span  across  lamina  boundaries.  Hence,  because  of 
the  laminated  character  of  the  material,  the  mesh  should  be  highly  refined  even 
where  the  strain  gradients  are  small. 


The  expense  of  modeling  each  lamina  individually  rapidly  becomes  intolerable 
as  the  number  of  laminae  increases.  Reference  1  presents  an  approximate  technique 
to  reduce  costs.  Laminate  theory  is  used  to  obtain  effective  extensional  moduli 
for  a  group  of  laminae.  Then  the  group  of  laminae,  rather  than  the  individual  lam¬ 
ina,  are  modeled  using  finite  elements.  This  approach  ignores  stacking  sequence 
effects  within  the  lamina  group.  Therefore,  the  flexural  and  flexural-extension 
coupling  properties  of  the  lamina  group  cannot  be  faithfully  modeled.  Reference  2 
presents  a  hybrid  analysis  for  thick  laminates.  In  this  analysis  (which  is  not 
finite  element  analysis)  the  laminate  is  divided  into  global  and  local  regions. 

The  terms  global  and  local  refer  to  the  detail  with  which  the 
modeled;  the  local  region  is  modeled  with  much  greater  detail  than  the  glob 
region.  Conceptually,  this  is  similar  to  using  a  finite  element  model 
smaller-  or  higher-order  elements  in  one  region  than  in  another.  However,  the 
Lalysis  in  re^rence  2  does  not  offer  the  inherent  flexibility  of  the  finite  ele¬ 
ment  method  for  modeling  complicated  geometries  and  for  performing  convergence 
checks.  The  objective  of  this  paper  is  to  introduce  a  new  type  of  two-dimensional 
(i.e.,  plane  stress  or  plane  strain)  finite  element  for  analysis  of  laminated  com- 

posites . 


The  element  is  a  four-node,  bilinear,  rectangular  element.  An  ordinary  bi¬ 
linear  rectangle  performs  poorly  in  modeling  bending-type  deformation.  The 
performance  can  be  improved  by  using  reduced  numerical  integration  or  substitute 
shape  functions  (ref.  3).  Because  of  the  multiple  laminae  within  the  element, 
numerical  integration  is  not  appropriate.  Therefore,  substitute  shape  functions 
are  used  to  improve  the  performance.  Explicit  integration  of  the  element  stiffness 
matrix  in  terms  of  generalized  displacements  minimizes  the  algebraic  effort 
required  to  account  for  the  various  laminae  within  a  single  element. 

After  describing  the  theoretical  aspects  of  the  element,  results  from  analyses 
of  several  simple  configurations  are  discussed. 
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SIMPLE  RECTANGULAR  ELEMENT  FOR  ANALYSIS  OF 
LAMINATED  COMPOSITE  MATERIALS 


John  D.  Whitcomb 
NASA  Langley  Research  Cehter 


ANALYSIS  OF  THICK  LAMINATES 


^  material  property  discontinuities 
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COMPARISON  OF  TRADITIONAL  AND  RLC  ELEMENT 


P 

P 


^Ex'lxlO^ 

^Ex*20xl0® 


P 


P 


X 


P 

P 


Traditional  element 


P 

P 


RLC  element 


RECTANGULAR  LAMINATED  COMPOSITE  ELEMENT 


X  Plane  stress/plane  strain 

X  Includes  stacking  sequence  effects 

X  Extensional  and  shear  stiffness 
X  Flexural  stiffness 
X  Extension-flexure  coupling 

^  Technique  valid  for  3-0  element 
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ELEMENT  CONFIGURATION 


X  Four  node  rectangle 
K  Lamina  interfaces  parallel  to  X-axis 


STIFFNESS  MATRIX 


c  . 


dA 


NORMAL  STRAINS: 

u«  a  +  bx  +  cy  +  dxy 
v«  a  +  bx  +  By  +  dxy 


€,-  -  b  +  dy 

I7  '  +  3x 


SHEAR  STRAIN: 

u»  e  +  fx  +  gy 
v»  e  +  fx  +  §y 


"xy 


9u  .  3v 


g  +  f  =  h 
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STIFFNESS  COEFFICIENTS 


CVI 

K 

13 

0 

0 

All® 

1  Cj,  (y 

1+1 

Vl) 

Kgg 

K 

23 

0 

0 

K 

33 

0 

0 

Aig  = 

I  C}g  (y 

i+i 

Vi) 

SYM 

0 

44 

K55 

Agg  = 

1  C^g  (y 

1+1 

Vl) 

A33= 

I  Cjg  (y. 

i+i 

Vi) 

II 

2 

lx 

tAii 

Bii= 

1/2Z  C}j 

(»!.. 

-y?) 

Ki2= 

*^13“ 

2 

2 

lx 

lx 

tBii 

tAig 

Bi2= 

I/2Z  Cjg 

(vf*. 

-vf) 

11 

OJ 

at 

2 

lx 

tDii 

Dii= 

1/31  C}, 

(Vt,. 

-  y?l 

^23~ 

2 

lx 

tB,2 

^33= 

2 

lx 

tAgg 

II 

2/3 

iJtAgg 

II 

to 

in 

2 

lx 

1-^33 

TRANSFORMATION  TO  NODAL  DISPLACEMENTS 


GENERALIZED  DISPLACEMENTS:  [K]  [  h]  =  IF] 

NDDAL  DISPLACEMENTS:  [K]  [  6]  =  [F] 

WHERE  [K]  =  [T]  ^  [K]  [T] 

[F]  =  [T]  ^  [F] 


TRANSFDRMATIDN  MATRIX: 


[A] 


BUT  [A]  = 


^NORMAL 


^SHEAR 


m  [fi] 


m= 


Jn_ORJ4)U. 

"^SHEAR 
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EXAMPLE  PROBLEMS 


MECHANICAL  LOAD 

X  THERMAL  LOAD 

^  ISOTROPIC 

X  DIFFERENT  MODULI  OR 

EXPANSION  COEFFICIENTS 


CANTILEVERED  BEAM  WITH  TIP  LOAD 

Laminate 
O  (H/H/H/H) 


Mesh 
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UNSYMMETRIC  CANTILEVERED  BEAM,  WITH 

THERMAL  LOAD 


€0—  nn  4 

€5=  .001 

€°=  -.001 
€°=  .001 


FINITE  ELEMENT  MODEL 


ERROR  IN  CALCULATED  TIP  DISPLACEMENT  <  .1% 


BEAM  UNDER  PURE  BENDING 


Y 


ASSUMED: 

€^=  b  +  dy 

€y=  c  +  dx  (CONSTANT) 
h  (ZERO) 


THEORETICAL: 

€jj=  b  +  dy 
€y=  -y(b  +  dy) 
€  =  ZERO 
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CANTILEVERED  BEAM  WITH  TIP  LOAD 


Laminate 
O  •  (H/H/H/H) 


Mesh 


SUMMARY 


^  Need  new  element  to  analyse  thick  laminates 


^  Developed  Rectangular  Laminated  Composite  element 

which  accounts  for  stacking  sequence  effects 

^  Performed  well  in  initial  tests 
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FORMULATION  OF  LAMINATED  BEAM  AND  PLATE  ELEMENTS  FOR  A  MICROCOMPUTER 

A.  Chen  and  T.  Yang 
Purdue  University 

Material  not  received  in  time  for  publication. 
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APPROXIMATE  ANALYSIS  OF  COMPOSITES  WITH  DAMAGE 
G.  Sendeckyj 

Flight  Dynamics  Laboratory 


Material  not  received  in  time  for  publication. 
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A  Method  of  Laminate  Design 


Ippei  Susuki 

Japan  National  Aerospace  Laboratory 
1 880  Jindaiji-machi,  Chofu,  Tokyo  1 82,  Japan 
presently 

Air  Force  Wright  Aeronautical  Laboratories/MLBM 
Wright-Patterson  Air  Force  Base.  Ohio  45433 


abstract 

The  objective  of  this  research  is  to  develop  a  design  method  to  translate 
the  strength  characteristics  of  unidirectional  composites  to  that  of  a  multi¬ 
directional  laminate  with  maximum  strength.  Composite  laminates  have  the 
advantage  that  their  mechanical  properties  can  be  tailored  through  fiber 
orientations,  ply  thickness  ratios,  and  stacking  sequences.  This  capability 
gives  a  designer  an  additional  degree  of  flexibility  to  obtain  the  desired 
stiffness  or  strength. 

This  work  is  limited  to  laminates  which  are  symmetric,  and  consist  of 
plies  with  0,90,45,  and  -45  degree  fiber  orientations.  A  method  is  developed 
to  find  the  optimal  combination  of  these  four  ply  groups  to  obtain  the 
maximum  strength  for  a  given  set  of  in-plane  loadings.  The  classical 
laminated  plate  theory  is  used  in  analyzing  the  stress  and  strain  of  each  ply. 
The  laminate  strength  is  determined  by  the  first  ply  failure  envelope  based  on 
the  quadratic  failure  criterion. 

Numerical  results  using  T300/5208  graphite/epoxy  composite  include: 

1 .  the  strength  of  optimized  laminate  is  about  six  times  as  high  as  that  of 
quasi-isotropic  one, 

2.  the  principal  stress  design  is  not  always  the  best, 

3.  theoretical, and  fractinal  ply  ratios  must  be  rounded  off  for  practical 
laminates.  The  effects  of  the  unit  ply  thickness  and  the  absolute  values  of 
applied  loads  will  be  illustrated. 

4.  the  penalty  of  using  balanced  laminates  under  non  principal  loads(with 
shears)  are  also  discussed. 

User  friendly  programs  on  PRIME(by  FORTRAN)  or  EPSON  QX-1 0(by  BASIC) 
for  this  work  are  available  to  interested  users  upon  request. 
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COORDINATE  SYSTEM  AND 
NOTATION  OF  MULTIDIRECTIONAL  LAMINATE 

S1,S2.S6:  LAMINATE  STRESSES 

(SHOWN  ARE  POSITIVE) 

<^k  :  PLY  ANGLE  OF  k-TH  PLY  GROUP 

v{</)k)  :  PLY  THICKNESS  RATIO  OF 

k-TH  PLY  GROUP 

H  :  ONE-HALF  OF  THE  TOTAL  THICKNESS 

hk  :  ONE-HALF  OF  k-TH  PLY  GROUP  THICKNESS 

Z 


ALUMINUM  AND  Gr/Ep  ELASTIC  PROPERTIES 


MATERIAL 

DENSITY 

Ikg/m^l 

Ex 

IGPal 

Ey 

IGPal 

Es 

IGPal 

Gr/Ep 

(T300/5208) 

1600 

181 

10.3 

0.28 

7.17 

ALUMINUM 

(2024-T3) 

2771 

12 

72 

0.30 

27.7 

STRENGTH  PROPERTIES 


LONG 

TENS 

LONG 

COMP 

TRANS 

TENS 

TRANS 

COMP 

SHEAR  IMPal 

Gr/Ep 

(T300/5208) 

1500 

1500 

40 

246 

68 

ALUMINUM 

(2024-T3) 

345 

345 

345 

345 

199 
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STRENGTH  SURFACE  OF 

A  QUASI-ISOTROPIC  LAMINATE  IN  REGION  T  {Sl=l  MPa) 


STRENGTH  SURFACE  OF 

A  QUASI-ISOTROPIC  LAMINATE  REGION  C  (Sl=-1  MPa) 
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LAMINATE  STRAINS  OF  (0p/90q)  LAMINATE 

(S1,S2,S6)  ■  (l,-.2,0)  MPa 
E1,E2:  NORMAL  STRAINS 
E6  :  SHEAR  STRAIN 

(IN  THIS  CASE  E6=0) 


^ - 1 - 1 - h-  v(0) 


R=  904.2  f 


STRENGTH  CURVE  OF  (0p/90q)  LAMINATE 
(S1.S2,S6)  -  (l,-.2,0)  MPa 


STRENGTH  SURFACE  OF 

(90q/45r/-45s)  LAMINATE  UNDER  (S1,S2,S6)=(1,.5,.3)  MPa 
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TETRAHEDRON,  SCHEMATIC  EXPRESSION  OF 
CONSTRAINT  ON  (0p/90q/45r/-45s)  LAMINATE  ANALYSIS 


SUBREGION 

PLY  THICKNESS  RATIO 

LAMINATE 

1 

v(0) 

fm 

v(45) 

v(-45) 

■ 

P 

Hi 

r 

S 

0 

A 

1 

0 

0 

0 

[0] 

1 

VERTEX 

B 

0 

1 

0 

0 

[90] 

2 

C 

0 

0 

1 

0 

[45] 

3 

0 

0 

0 

0 

1 

[-45] 

4 

AB 

l-x 

X 

0 

0 

[0p/90q] 

5 

AC 

l-x 

0 

X 

0 

[0p/45r] 

6 

EDGE 

OA 

l-x 

0 

0 

X 

[0p/-45s] 

7 

BC 

0 

l-x 

X 

0 

8 

OB 

0 

■BB 

0 

X 

E 

OC 

0 

0 

l-x 

X 

[45  r/- 45s] 

m 

OBC 

0 

HDBI 

X 

y 

[90q/45r/-45s] 

n 

SURFACE 

OPC 

IDSi 

0 

X 

y 

[Op/45r/-45s] 

m 

OAB 

X 

0 

y 

[0p/90q/-45s] 

E 

ABC 

HB3 

X 

y 

0 

[0p/90q/45r] 

IS 

BODY 

0-ABC 

X 

y 

Z 

(0p/90q/45r/-45s] 

0  <  X,  y,  z  <  I 
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STRENGTH  SURFACE  OF 

OPTIMIZED  LAMINATE  IN  REGION  T  (S1=1  MPa) 


PLY  GROUPS  OF  OPTIMIZED  LAMINATE  IN  REGION  T  (Sl=l  MPa) 
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STRENGTH  VARIATION  IN 
OPTIMIZED  LAMINATE  AT  S6=.l  MPa  IN  REGION  T 


PLY  THICKNESS  RATIOS  OF  OPTIMIZED  LAMINATE  AT  S6=.l  MPa  IN  T 
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C  13  Stress  resultant  N1  N2  N6  (fIN/tn)  i  -1«000  0.000  0.000 

flaterxal  {T300/5208  Graphite/Epoxy  with  f*-0.5 
No.  NTPY  H(fflm)  NHPY  C  0  90  45-453  FPTF  Rmln  R(  0)R(  90)R(  45)R(-45) 


1 

6 

0,751 

3 

( 

3 

0 

0 

0) 

0 

1.1251 

1.125 

0 . 000 

0.000 

0 . 000 

2 

6 

0.75! 

3 

( 

2 

1 

0 

0) 

0 

1.060! 

1.060 

2,186 

0 . 000 

0 . 000 

3 

8 

1.001 

4 

( 

3 

1 

0 

0) 

0 

1.552! 

1.552 

3.220 

0.000 

0.000 

4 

3 

1.00! 

4 

( 

4 

0 

0 

0) 

0 

1.5001 

1.500 

0.000 

0.000 

0.000 

5 

8 

1.00! 

4 

( 

2 

2 

0 

0) 

0 

1.108! 

1.108 

2.269 

0.000 

0,000 

6 

10 

1.25! 

5 

( 

4 

1 

0 

0) 

0 

2.031! 

2.031 

4.238 

0.000 

0 . 000 

7 

10 

1.25! 

5 

( 

5 

0 

0 

0) 

0 

1.875! 

1.875 

0.000 

0.000 

0.000 

8 

10 

1.25! 

5 

( 

3 

2 

0 

0) 

0 

1.618! 

1.618 

3.326 

0,000 

0,000 

9 

10 

1.25! 

5 

( 

3 

1 

0 

1) 

0 

1.339! 

1.339 

2.680 

0.000 

1.871 

10 

10 

1.25! 

5 

( 

3 

1 

1 

0) 

0 

1.339! 

1.339 

2.680 

1.871 

0.000 

nextiSl 


t  Nl(nN/m)«-10 

Uni-axial  Compression 

C  23  Stress  resultant  N1  N2  N6  (PIN/m)  :  -10.000  0.000  0.000 

Material  tT300/5208  Graphite/Epoxy  with  f*-0.5 
No.  NTPY  H(mffl)  NHPY  C  0  90  45-453  FPTF  Rmin  R(  0)R(  90)R(  45)R(-45) 


1 

48 

6.00! 

24 

( 

22 

2 

0 

0) 

0 

1.0321 

1.032 

2.210 

0 . 000 

0.000 

2 

48 

6.00! 

24 

( 

21 

3 

0 

0) 

0 

1.023! 

1.023 

2.164 

0,000 

0.000 

3 

48 

6.00! 

24 

( 

23 

1 

0 

0) 

0 

1.008! 

1.008 

2.207 

0.000 

0.000 

4 

50 

6.25! 

25 

( 

23 

2 

0 

0) 

0 

1.074! 

1.074 

2.304 

0.000 

0 . 000 

5 

50 

6.25! 

25 

( 

22 

3 

0 

0) 

0 

1.068! 

1.068 

2.262 

0.000 

0,000 

6 

50 

6.25! 

25 

( 

24 

1 

0 

0) 

0 

1.048! 

1.048 

2.296 

0,000 

0.000 

7 

50 

6.25! 

25 

( 

21 

4 

0 

0) 

0 

1.046! 

1.015! 

1.046 

2.196 

0.000 

0.000 

8 

50 

6.25! 

25 

( 

20 

5 

0 

0) 

0 

1.015 

2.119 

0 . 000 

0 . 000 

9 

52 

6.50! 

26 

( 

24 

2 

0 

0) 

0 

1.1171 

1.117 

2.398 

0.000 

0 . 000 

10 

52 

6.50! 

26 

( 

23 

3 

0 

0) 

0 

1.113! 

1.113 

2.359 

0.000 

0.000 

next? 

54 

6.75 

27 

( 

27 

0 

0 

0) 

0 

1.013 

C 


43  Stress  resultant  N1  N2  N6  (MN/m)  :  2.000  0.400 


0.000 


Material  tT300/5208  Graphite/Epoxy  with  f=-0.5 
No.  NTPY  H(mm)  NHPY  C  0  90  45-453  FPTF  Rmin  R(  0)R(  90)R(  45)R(-45) 


1 

28 

3.50! 

14 

( 

8 

0 

3 

3) 

45 

1.025! 

1.339 

0.000 

1.025 

1.025 

2 

30 

3.75! 

15 

( 

9 

0 

3 

3) 

45 

1.101! 

1.382 

0.000 

1.101 

1.101 

3 

30 

3.75! 

15 

( 

8 

0 

3 

4) 

-45 

1.039! 

1.465 

0.000 

1.111 

1.039 

4 

30 

3.75! 

15 

( 

8 

0 

4 

3) 

45 

1.039! 

1.465 

0.000 

1.039 

1.111 

5 

30 

3,75! 

15 

( 

7 

0 

4 

4) 

45 

1.024! 

1.523 

0 . 000 

1.024 

1.024 

6 

30 

3.75! 

15 

( 

10 

0 

2 

3) 

-45 

1.005! 

1.191 

0.000 

1.115 

1.005 

7 

30 

3.75! 

15 

( 

10 

0 

3 

2) 

45 

1.005! 

1.191 

0.000 

1.005 

1.115 

8 

30 

3.75! 

15 

( 

11 

0 

2 

2) 

-45 

1.002! 

1.045 

0.000 

1.002 

1.002 

9 

32 

4.00 ! 

16 

( 

10 

0 

3 

3) 

45 

1.171! 

1.417 

0.000 

1.171 

1.171 

10 

32 

4.00! 

16 

( 

8 

0 

4 

4) 

45 

1.127! 

1.615 

0.000 

1.127 

1.127 

11 

32 

4.001 

16 

( 

9 

0 

3 

4) 

-45 

1.121! 

1.522 

0.000 

1.201 

1.121 

12 

32 

4.00! 

16 

( 

9 

0 

4 

3) 

45 

1.121! 

1.522 

0.000 

1.121 

1.201 

13 

32 

4.00! 

16 

( 

11 

0 

3 

2) 

45 

1.054! 

1.215 

0.000 

1.054 

1.169 

14 

32 

4.001 

16 

( 

11 

0 

2 

3) 

-45 

1.054! 

1.215 

0.000 

1.169 

1.054 

15 

32 

4.00! 

16 

( 

8 

0 

3 

5) 

-45 

1.047! 

1.546 

0.000 

1.171 

1.047 

16 

32 

4.00! 

16 

( 

8 

0 

5 

3) 

45 

1.047! 

1.546 

0.000 

1.047 

1.171 

17 

32 

4.00! 

16 

( 

12 

0 

2 

2) 

45 

1.042! 

1.065 

0.000 

1.042 

1.042 

18 

32 

4.00! 

16 

( 

7 

0 

4 

5) 

-45 

1.031! 

1.606 

0.000 

1.077 

1.031 

19 

32 

4.00! 

16 

( 

7 

0 

5 

4) 

45 

1.031! 

1.606 

0.000 

1.031 

1.077 

20 

32 

4.00! 

16 

( 

10 

0 

2 

4) 

-45 

1.028! 

1.301 

0.000 

1.221 

1.028 

21 

32 

4.00! 

16 

( 

10 

0 

4 

2) 

45 

1.028! 

1.301 

0.000 

1.028 

1.221 

22 

32 

4.001 

16 

( 

11 

1 

2 

2) 

90 

1.028! 

1.315 

1.028 

1.143 

1.143 

23 

32 

4.00! 

16 

( 

12 

2 

1 

1) 

90 

1.027! 

1.253 

1.027 

1.122 

1.122 

24 

32 

4.00! 

16 

( 

13 

3 

0 

0) 

90 

1.026! 

1.211 

1.026 

0.000 

0.000 

25 

32 

4.001 

16 

( 

13 

1 

1 

1) 

90 

1.0171 

1.051 

1.017 

1.033 

1.033 

26 

32 

4.00! 

16 

( 

14 

2 

0 

0) 

90 

1.015! 

1.042 

1.015 

0.000 

0.000 
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C  13  Stress  resultant  N1 
Material  :T300/5208 
No.  NTPY  H(niffl)  NHPY  C 

1  10  1.251  5  ( 

2  12  1.50i  6  ( 

3  12  1.501  6  ( 

4  14  1.75!  7  ( 

5  14  1.751  7  ( 

6  14  1.75!  7  ( 

7  16  2.00!  8  ( 

8  16  2.00!  8  ( 

9  16  2.00!  8  ( 

10  16  2.00!  8  ( 

C  23  Stress  resultant  N1 
Material  JT300/5208 
Mo.  NTPY  NHPY  C 


N2  N6  (MN/m)  i  1.000 
Graphite/Epoxy  with  f" 
0  90  45-453  FPTF  Rmin 


4 

0 

1 

0) 

45 

1.033! 

5 

0 

1 

0) 

45 

1.270! 

4 

0 

2 

0) 

45 

1.032! 

6 

0 

1 

0) 

45 

1.490! 

5 

0 

2 

0) 

45 

1.297! 

4 

0 

3 

0) 

45 

1.020! 

7 

0 

1 

0) 

0 

1.605! 

6 

0 

2 

0) 

45 

1.560! 

5 

0 

3 

0) 

45 

1.283! 

5 

0 
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SUPPRESSION  OF  DELAMINATION  IN  COMPOSITE  LAMINATES  SUBJECTED  TO  IMPACT  LOADING 

C.  T.  Sun 

School  of  Aeronautics  and  Astronautics 

Purdue  University 

West  Lafayette,  Indiana  47907 


ABSTRACT 


Two  major  failure  modes  are  predominant  in  composite  laminates  subjected  to  low 
velocity  impact,  i.e.,  matrix  cracks  and  delamination  [1].  Delamination,  in  particular, 
may  cause  a  great  reduction  in  compressive  strength  and,  thus,  suppression  of  delamina- 
tion  is  highly  desirable. 

Analysis  results  indicate  that  delamination  is  induced  by  matrix  cracks.  Two  types 
of  matrix  cracks  are  induced  by  transverse  shear  stresses  and  bending  stresses.  Bending 
cracks  usually  occur  in  the  bottom  layer  where  bending  stresses  are  the  greatest,  while 
transverse  shear  cracks  occur  more  often  in  the  interior  layers.  Delamination  can  re¬ 
sult  from  branching  of  both  bending  cracks  and  transverse  shear  cracks.  In  thin  lami¬ 
nates  it  is  not  produced  by  reflection  of  waves  propagating  through  the  thickness  of  the 
1  ami na te . 

Two  methods  of  delamination  suppression  are  investigated:  adding  through-the- thi ck- 
ness  reinforcements  by  stitching,  and  adding  tough  adhesive  layers  along  interfaces  of 
the  laminate.  Experimental  results  indicate  that  stitching  can  reduce  the  spreading  of 
delamination  cracks  but  cannot  prevent  initiation  of  delamination.  On  the  otherhand, 
the  use  of  adhesive  layers  proves  to  be  quite  effective  in  suppressing  del  ami nati on . 

Thin  structural  adhesives  were  placed  along  the  interfaces  of  a  composite  laminate 
subjected  to  impact  loading.  The  base  line  specimen  was  AS4/3501-6  graphite/epoxy 
[Or/OOr/Or]  laminate.  One  group  of  specimens  contained  layers  of  5  mil  thick  adhesive 
film  (NlOOO  by  American  Cyanamid)  between  the  0°-  and  90°-plies.  Impacted  beam  speci¬ 
mens  were  sectioned  transversely  and  longitudinally  through  the  impact  site  and  were 
examined  using  an  optical  electron  microscope  and  then  photographed.  From  these  photo¬ 
graphs,  matrix  cracks  and  delamination  size  were  measured.  The  X-ray  technique  was  also 
used  to  map  the  delamination  area. 

When  compared,  the  base  line  laminate  and  the  one  with  adhesive  layers  exhibited 
significant  differences  in  damage  mode.  For  instance,  the  threshold  velocity  at  which 
damage  (matrix  cracking  or  delamination)  occurred,  was  found  to  be  much  lower  for  the 
laminate  with  no  adhesives.  It  was  found  that  a  delamination  crack  could  not  be  initia¬ 
ted  in  the  laminate  with  adhesives  until  the  impact  velocity  reached  26  m/sec,  while 
under  a  13  m/sec  impact,  the  laminate  without  adhesives  began  to  suffer  del  ami nati on . 

Transverse  shear  crack  in  the  upper  lamina  was  also  noted  to  be  affected  by  the 
presence  of  adhesive  layers.  For  example,  at  15  m/sec  transverse  shear  cracks  appeared 
in  the  upper  lamina  in  the  specimen  without  adhesives,  while  no  transverse  cracks  emerged 
in  the  upper  lamina  in  the  beam  with  adhesives  at  velocities  below  26  m/sec.  There  was 
indication  that  delamination  could  be  caused  by  transverse  cracks. 

The  adhesive  layers  in  the  laminate  appeared  to  have  affected  the  size  of  indenta¬ 
tion  during  impact.  As  a  result,  the  distribution  of  contact  pressure  was  also  affected. 
Experiments  were  conducted,  and  the  relations  between  the  contact  area  and  impact  veloci¬ 
ty  for  laminates  with  adhesives  and  without  adhesives  were  obtained.  The  results  showed 
that  adhesive  layers  tended  to  enlarge  the  contact  area  and,  as  a  result,  spread  out  the 
contact  pressure.  The  distribution  of  the  contact  pressure  was  found  to  significantly 
affect  the  transverse  shear  stress  distribution  through  the  thickness.  This  result  was 
used  successfully  to  interpret  the  difference  in  transverse  shear  crack  patterns  in  the 
two  types  of  laminates  due  to  impact  loading. 

The  effect  of  adhesive  layers  on  transverse  shear  stress  distribution  was  studied 
by  using  2-D  finite  elements.  The  results  show  that  the  adhesive  layer  can  diffuse  shear 
stress  concentration  in  the  top  layer  (impacted  side)  resulting  in  a  significant  reduc¬ 
tion  in  transverse  shear  cracks.  Effects  of  adhesive  layer  locations  were  also  invest!- 
gated . 
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OBJECTIVE 

TO  INVESTIGATE  THE  EFFECTIVENESS  IN 
USING  STITCHING  AND  ADHESIVE  LAYERS 
FOR  SUPPRESSION  OF  DELAMINATION 


APPROACH 

»  IMPACT  [O5/9O5/O5]  LAMINATES  WITH 
STITCHING  AND  ADHESIVE  LAYERS 

"  DETERMINE  CONTACT  AREAS 

“  USE  PHOTOMICROGRAPHS  TO  OBSERVE 
IMPACT  DAMAGE 

»  USE  2-D  FINITE  ELEMENTS  TO  PERFORM 
DYNAMIC  STRESS  ANALYSIS  TO  INTERPRET 
EXPERIMENTAL  RESULTS 


crack 


TRANSVERSE  SECTION 
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(b)  [O5/A/9O5/A/O5] 

Transverse  cross  sections  of  laminates  after  inpact 
of  a  1/2"  steel  ball  at  velocity  21  tn/s 


Transverse  cross  sections  of  laminates  after  impact 
of  a  1/2"  steel  ball  at  velocity  about  26  m/s 


0 

0 


UPPER  INTERFACE  OELAMINATION 
LOWER  interface  DELAMINATION 


5  MM 

IMPACT  velocity  16m/ sec 

Delaminaiion  zones 


UPPER  INTERFACE  DELAMINATION 
LOWER  INTERFACE  OELAMINATION 
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a  =  0.10  mm 
b  =  0.i|5  mm 
c  =  0.67  mm 
d  =  1.25  nm 


u loi  . uu I 

[S05/A/O5/A/S0S] 
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CONTACT  LENGTH  (nm) 


impact  velocity  34  m/sec  impact  velocity  14.6  m/sec 
delamination  in  05-layer  adhesive  layers  in  O5  layer 
fiber  breakage 


impact  velocity  14  m/sec 
adhesive  layer  near  top 
and  bottom  surface 


Impact  velocity  14.6  m/sec 
adhesive  layer  at  the 
middle  surface 
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CONCLUSIONS 

“  STITCHING  CAN  REDUCE  DEL/V1INATION  SIZE 
BUT  CANNOT  SUPPRESS  INITIATION. 

^  ADHESIVE  LAYERS  WHEN  PLACED  AT  INTERFACES 
ARE  EFFECTIVE  IN  SUPPRESSING  DELAMINATION. 

"  ADHESIVE  LAYERS  CAN  REDUCE  TRANSVERSE 
SHEAR  CRACKS. 

°  ADHESIVE  LAYERS  CANNOT  STOP  BENDING  CRACKS. 
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Gwynedd  Plaza  II 
Bethlehem  Pike 

Spring  House,  Pennsylvania  19477 


ABSTRACT 


Structural  materials  in  space  are  subjected  to  severe  thermal  cycling  as  the 
structure  moves  between  full  sunlight  and  the  shadow  of  the  earth.  The  range  of 
temperatures  experienced  can  be  as  large  as  1160k  to  5890k.  For  polymeric  matrix  fiber 
composites  this  produces  creep  deformations  and  relaxation  of 

accurately  assessed.  The  purpose  of  this  study  is  to  (a)  provide  analytical  methods  t 
determine  the  thermoviscoelastic  (TVE)  properties  of  graphite/polymer  unidirectional 
fiber  composite  (UFC),  (b)  to  apply  this  information  to  the  assessment  of  composite 
Structural  response. 

TVE  characterization  of  polymeric  matrix  materials  is  not  a  simple  matter  and 
various  schemes  have  been  devised  in  addressing  this  problem.  In  this  study  we  have 
adopted  the  simplest  representation  which  can  adequately  characterize  the  response  o  a 
TVE  polymer.  The  representation  used  retains  the  temperature  dependence  of  the  initial 
elastic  modulus  while  the  time  dependent  response  is  governed  by  a  temperature  dependent 
time  shift  function. 

It  will  be  recalled  that  for  isothermal  viscoelasticity  there  is  a  correspondence 
principle  which  directly  relates  effective  viscoelastic  and  elastic  properties  for  any 
model  of  a  composite  materials  [1,2]  but  no  such  correspondence  exists  for  a  composite 
with  a  TVE  constituent.  It  therefore  becomes  necessary  to  resort  to  direct  analysis  ot 
a  model.  We  have  employed  two  models  of  a  UFC  -  the  composite  cylinder  assemblage  (CCA) 
model  which  is  suitable  for  axisymmetric  states,  including  free  thermal  expansion  and  for 
axial  shear,  and  the  hexagonal  array  model  for  the  case  of  transverse  shear.  It  may  be 
shown  that  any  average  stress  or  strain  state  can  be  decomposed  into  axisymmetric  part, 
transverse  shear  and  axial  shear.  Thus  analysis  of  these  three  states  provides  the 
general  stress-strain  response  of  a  UFC. 

The  advantage  of  the  CCA  model  is  in  that  the  analysis  reduces  to  that  of  a  single 
composite  cylinder.  Carrying  out  such  analyses  for  axisymmetric  states  and  axial  shear 
it  was  found  that  in  each  case  the  problem  was  reduced  to  solution  a  system  of  integral 
eauations  in  the  time  which  is  easily  done  numerically  in  an  incremental  fashion.  The 
CCA  model  is  not  amenable  to  transverse  shear  analysis  and  therefore  this  case  has  been 
analyzed  in  terms  of  the  hexagonal  array  model  with  TVE  finite  elements  which  were 
specially  developed  for  this  study.  Note  that  the  CCA  and  the  hexagonal  array  are  both 
transversely  isotropic  models  as  required  and  that  their  predictions  for  elas  ic 
constants  of  UFC  are  numerically  literally  identical. 


The  constitutive  relations  for  UFC  developed  have  to  be  utilized  to  study  the 
vibrations  of  a  cantilever  beam.  With  respect  to  this  vibration  problem  it  is  important 
to  recognize  that  a  UFC  will  have  small  time  dependent  response  when  loaded  in  fiber 
direction  and  large  time  dependence  when  sheared  along  the  fiber  direction  (axial  shear) . 
Thus,  in  bending  vibrations  of  a  beam  which  is  uniaxially  reinforced  in  the  direction  of 
the  beam  axis  the  effect  of  shear  is  very  significant  and  therefore,  we  have  based  the 
analysis  on  the  Timoshenko  beam  concept  which  takes  into  account  shear  deformation  and 
rotatory  inertia.  For  this  purpose  the  Timoshenko  beam  vibration  equations  have  been 
converted  to  viscoelasticity  in  terms  of  complex  moduli  (axial  Young's  modulus  and  axial 
shear  modulus)  of  the  UFC.  These  have  been  evaluated  at  different  temperatures  in  terms 
of  matrix  complex  moduli  and  fiber  elastic  moduli.  We  have  evaluated  the  attenuation  of 
free  bending  vibrations  of  a  cantilever  with  circular  tubular  section  without  and  with 
the  shear  deformation  effect.  The  results  show  that  shear  deformation  provides  very 
significant  attenuation  of  vibrations.  Indeed  it  is  seen  that  modes  higher  than  the  first 
are  effectively  eliminated  by  the  damping. 
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OBJECTIVES 


o  OBTAIN  THERMO-VISCOELASTIC  CONSTITUTIVE  RELATIONS 
FOR  UNIDIRECTIONAL  FIBER  COMPOSITES 

o  INVESTIGATE  COMPOSITE  STRUCTURAL  RESPONSE 


APPROACH 


O  MICROMECHANICS  MODELLING 


o  STRUCTURAL  MODELLING 


THERMO-VISCOELASTIC  RESPONSE  OF  UNIDIRECTIONAL  FIBER  COMPOSITE 


FIBERS  -  THERMO-ELASTIC  ANISOTROPIC  GRAPHITE,  CARBON 

MATRIX  -  THERMO-VISCOELASTIC  POLYMER 
INPUT  -  (Jilt)  find  -  ^ij  (t) 

RELAXATION 

INPUT  -  (|>(t)  find  -  e^j(t) 

CREEP 
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THERMO-ELASTIC  DILATATION  RESPONSE 


a. .  =  06 . .  +  s. . 
1:  1]  ID 


0  5 


=  e6..  +  e.. 
ID  ID  ID 


^  3  ^^ll‘*'^22'''^33’ 


0  [(|>(t)]  =  3K  [((.(t)]  |e(t)  -  a[<l)(t)]  ().{t) 

THERMO-ELASTIC  RELATION 

<|)(t)  -  TIME  DEPENDENT  TEMPERATURE 


THERMO-VISCOELASTICITY 


ISOTHERMAL  ifi  =  CONST 

rt  ae.. 

S^j(t)  =  2G(t,<j>)  e^j(O)  +  2y^+  G(t-t’,((.)  -^df 

eij(t)  =|g(t,<t.)  s.j(o)  g(t-t',<^)  <af 

G(t,(j))  -  SHEAR  RELAXATION  MODULUS 
g(t,(f)  -  SHEAR  CREEP  COMPLIANCE 

THERMORHEOLOGICALLY  SIMPLE 

rt 

S^j(t)  =  2G(^)  e^j(O)  +  2j^l  df 

as 

e.j(t)  =^(?)  S..(0)  g(C-C')  -5^  df 

reduced  time 

a  -  HORIZONTAL  SHIFT  FACTOR 

G(t,(|)^)  =  F(logt) 

G(t,(ti)  =F[log(t/a)] 
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THERMORHEOLOGICALLY  COMPLEX 


TRANSVERSELY  ISOTROPIC  FIBER  COMPOSITE 


COMPOSITE  CYLINDER  HEXAGONAL  COMPOSITE  CYLINDER 

ASSEMBLAGE  ARRAY  ASSEMBLAGE 
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AXISYmETRI 


ui  « 

u^{b,t)  -  e^{t)b 

“  2  ^^22  ^33^ 

“  2  (*^22  ^33) 


(|)  =  (t»(t) 


FREE  THERMAL  EX 


Matrix 


AXIAL  SHEAR 


PERIODIC  HEXAGONAL  ARRAY  AND  REPEATING  ELEMENT 
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ELASTIC  BEAM  VIBRATIONS 


NEGLECT  SHEAR  DEFORMATION  AND  ROTATORY  INERTIA 


WITH  SHEAR 

DEFORMATION  AND 

ROTATORY 

INERTIA. 

TIMOSHENKO 

BEAM 

El 

3x 
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TV  9  W 

+  pA  „ 

-  pi 

3t^ 

V  k'G 

W 


2  2 


-  SHEAR  COEFFICIENT 

VISCOELASTIC  VIBRATIONS 


W (x, t)  =  W (x) 

£•>£*(01))  +i  E"(co)  G  G*(aj)  +  i  G"  (w) 

1 

i  COMPLEX  NATURAL  FREQUENCY 


W^(x,t)  =  W^(x) 


Attenuation 


Mode 

Sinusoidal 
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FIRST  MODE  NATURAL  FREQUENCIES.  CANTILEVER  BEAM 


400  h 


w  300 


200 


100 


Bernoulli  -Euler  Beam 
Timoshenko  Beam 


L  -  40  in 
R  =  0.5  in 
t  =  0.01  in 


CARBON-EPOXY 


J- 


J_ 


_L 


-300  -200 


-100  0  100 
TEMPERATURE  {°F) 


200 


10"- 


10 


,10 


10 


300  400 


ENVELOPE  OF  FIRST  MODE  RESPONSE.  CANTILEVER  BEAM 


,22 


(SEC 


AMPLITUDE 


ENVELOPE  OF  SECOND  MODE  RESPONSE.  CANTILEVER  BEAM 


-  Bernoulli-Euler  Beam 

-  Timoshenko  Beam 


<\>  =  86“F 


CARBON-EPOXY 


.^00  200  300  400 


-0.5  -I 

I 

I 


TIME  (SEC) 


J  0.0 


SUMMARY 


o  TEMPERATURE  DEPENDENCE  OF  INITIAL  PROPERTIES  SIGNIFICANT 


o  SECONDARY  STRESS  DOMINATE  BEAM  DAMPING 


o  STRUCTURAL  MODELS  MUST  INCORPORATE  SECONDARY  STRESSES 
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COMPOSITE  MECHANICS/RELATED  ACTIVITIES  AT  LEWIS  RESEARCH  CENTER 

C.  C.  CHAMIS,  C.  A.  GINTY,  AND  P.  L.  N.  MURTHY 

NASA  LEWIS  RESEARCH  CENTER 
CLEVELAND,  OHIO  44135 


ABSTRACT 


Lewis  research  activities  and  progress  in  composite  mechanics  and  closely  related 
areas  are  summarized.  The  research  activities  summarized  include:  UJ  Composite 
Mechanics;  (2)  Computer  Programs  for  Composites;  and  (3)  High  Temperature  Composites. 
The  research  activities  focus  on:  (1)  Composite  mechanics  --  simplified  procedures 
for  sizing  laminates,  substructuring  in  composite  mechanics  for  free  edge  stress- 
field  behavior,  interlaminar  fracture  toughness  in  composites,  and  progressive  tr£ 


vsis,  therinoviscopiasT:ic  compubitc  ucuavxv^i,  i.  1 

blades  and  structural  tailoring  of  advanced  turboprops;  and  (3)  High  temperature 
composites  --  test  methods  development  and  characterization,  composite  burner  liner 
and  refractory  wire  reinforced  superalloys  for  SSME  turbopump  blades. 


0B3ECTIVE 


COMPOSITE  MECHANICS /RELATED  ACTIVITIES  SUMMARY  OF  LBJIS  RESEARCH  ACTIVniES  AM)  PROGRESS  IN: 

AT  LEWIS  RESEARCH  CENTER 

0  COMPOSIIE  MECHANICS 

0  COMPUTER  PROGRA/'IS  FOR  COMPOSITES 

0  HIGH  TEMPERATURE  COMPOSITES 

C.  C.  CHAMIS.  C.  A.  GINTY,  AND  P.  L.  N.  MURTHY 
NASA  LEWIS  RESEARCH  CENTER 
CLEVELAND.  OHIO  MH135 


COMPOSITE  MECHANICS 


ELEVENTH  ANNUAL  MECHANICS  OF  COMPOSITES  REVIEW 
DAYTON.  OHIO,  OCTOBER  22-24,  1985 


0  SIMPLIFIED  PROCEDURES  FOR  SIZING  LAMINATES 

0  SUBSTRUCTURING  IN  COMPOSITE  MECHANICS:  WIDTH.  LOADING  CONDITION  AND 
ENVTROfllENTAL  EFFECTS  ON  FREE  EDGE  STRESSES 

0  INTERLAMINAR  FRACTURE  IN  COMPOSITES 

0  PROGRESSIVE  FRACTURE  OF  COMPOSITES  WCTH  AND  WHHOUT  DEFECTS 
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SCHEMATICS  OF  SELECT  COMPOSITE  STRUCTURAL  COMPONENTS  WITH  RESPECTIVE 
GEOMETRY  AND  TYPICAL  LOADINGS 


C.  BOX  BEAM 


V-2452 


LEWIS  DEVELOPED  DESIGN  PROCEDURE 

CONSISTS  OF: 

•TABULATED  TYPICAL  PROPERTIES  OF  VARIOUS  UNIDIRECTIONAL 
COMPOSITES 

•  SUITABLE  GRAPHS  TO  EXPEDITE  REPETITIOUS  CALCULATIONS 

•  PLY  STRESS  INFLUENCE  COEFFICIENTS 

•  EXPLICIT,  SIMPLf  EQUATIONS  FOR  HYGROTHERMAL  EFFECTS  AND 

FATIGUE 

•APPROXIMATE,  SIMPLf  EQUATIONS  FOR  BUCKLING 


ELASTIC  PROPERTIES  OF  AS-GRAPHITE-FIBER/EPOXY  (AS/E)  ±0  LAMINATES 


REDUCED  STIFFNESSES  OF  AS  GRAPHITE-FIBER/EPOXY 
(AS/E)  ±e  LAMINATES 


POISSON'S  MODULI 

RATIO  (nipsi) 
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TliCRMAL  EXPANSION  COEl TICIENTS  OF  AS  GRAPlinE-FlBER/EPOXY 
(AS/E)  ±9  LAMINATES 


MOISTURE  EXPANSION  COEFFICIENTS 
OF  AS  GRAPIIITE-FiBER/EPOXY 

(AS/E)  ±fl  LAMINATES 


. lSxlO"2 


PLY  STRESS  INFLUENCE  COEFFICIENTS  (PSIC)  -  DEFINITION 

PSIC  Wl  RELATE  THE  8  PLY  MATERIAL  AXES  STRESSES  TO  LAMINATE  STRUCTURAL  AXES 
STRESSES 


LAMINATE  STRUCTURAL  AXES  (X,  Y)  STRESSES  9  PLY  MATERIAL  AXES  (1, 2)  STRESSES 

°je22  '  '^T/X  °cxx  *'^T/Y  ®cyy  '^•^T/S  °cxy 
®^12  '  '^S/X  °cxx  *  '^S/Y  °cyy  *  '^S/S  °cxy 


SAMPLE  DESIGN;  COMPOSITE  PANEL  SUBJECTED  TO  IN-PLANE  COMBINED  LOADS 


COMPOSITE  SYSTEM:  AS/E,  ABOUT  0.6  FVR 
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GENERAL  LAMINATE  DESIGN  RESULTS  SUMMARY 
[LAMINATE  CONFIGURATION:  [+45/0/ 90/0]  25  ^SIE;  t,.”  0.10  In.] 


MARGINS  OF  SAFETY  (M.  0.  S.  I  FOR: 


DISPLACEMENT 

PLY  STRESSES 

BUCKLING 

TYPE 

M.  0.  S. 

PLY 

M.  0.  S. 

C 

ASE,  ks 

M.  0.  S. 

u/a 

v/b 

A9 

0.43 

1.94 

.33 

0 

+45 

-45 

90 

2.77 

0.79 

4.43 

6.00 

'’/22 

0.61 

00 

^27 

.12 

Hu 

1.30 

00 

00 

1.30 

0 

40 

Ocyy 

0 

20 

^cxy 

20 

20 

-0.74 

3.35 

^AT  SPECIFIED  LOAD  (AT  DESIGN  LOAD  M.  S.  =  0. 381. 

CD-85-16564 


FREE-EDGE  SUPERELEMENT 

(224  F.E.  •.  2355  D.O.F.I 


CASES  STUDIED 

•  FREE-EDGE  AXIAL  STRESS  VARIATION  VERSUS  PLY  ANGLE 

•  WITII/THICKNESS  RATIO  EFFECTS  ON  FREE-EDGE  STRESSES 

•AXIAL  TENSION  -  REFERENCE 

•IN-PLANE  SIIEAR/BENDING 

•  OUI-OF-PLANE  SHEARIBENDING 

•  BENDING  MOMENT 

•TWISTING  MOMENT 

•COMBINED  LOADING 

(AXIAL  TENSION/IN-PLANE  SHEAR) 

•  UNIFORM  THERMAL  LOAD 
•UNIFORM  MOISTURE  LOAD 

CD-85-15810 


X 


FINITE  ELEMENT  MODEL  OF  ANGLEPLIED  LAMINATE 
UNDER  UNIFORM  EXTENSIONAL  STRESS 


Y 


- - 28" - - 

11 

y' 

y' 

y' 

X 

S 

■'-FREE  EDGE  SUPERELEMENT 


F.  E..  STATISTICS 
1589  20-NODE  SOLID  F.  E. 

224  F.  E.  's  IN  SUPERELEMENT 
22683  D.  0.  F. 
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PROPERTIES 

(AS/EP  0°  PLY  AND  INTERPLY  LAYER) 


Oil’ 

^11  ^12  *^13  ^14  ®15  ^16 

^11 

022 

G22  G23  G24  G25  G26 

£22 

033 

G33  G34  G35  G36 

X 

£33 

044 

G44  G45  G46 

£44 

055 

*^55  ^56 

^55 

”66 

^66 

^66 

^11 

19. 4  Gi2  '  0. 5164  G13  =  0. 5164  G22 

=  1.384 

G23 

0.4473  G33  =  1.483  G44  “  0.6667  G55 

=  0,3974 

G56  »  0,6667  E  -  0.50;  V  '  0.35 
(UNITS  X  10^  psi) 


COMPARISON  OF  FREE-EDGE  STRESSES 
DUE  TO  ENVIRONMENTAL  LOADS 

{±10®AS;EAPLAT0.55FVR) 

INTERLAMINAR 


COMPARISON  OF  FREE-EDGE  STRESSES  DUE  TO  MECHANICAL  LOADS 

(tlO^ASIE  APLAT0.55FVRI 


NORMALIZED 

STRESS, 

®2Z'®XZ*®XX.CLT 


WIDTH  TO  THICKNESS  RATIO  EFFECTS  OF 
INTERLAMINAR  FREE-EDGE  STRESSES 

(±45®AS/EAPL.  0.55FVRI 


128 


SCHEMATIC  OF  FLEXURAL  TEST  FOR  INIERLAMINAR 
FRACTURE  MODE  TOUGHNESS 


END  NOTCH-FLEXURE  ENERGY  RELEASE  RATECOMPARISONS 


GENERAL  PROCEDURE  FOR  PREDICTING  COMPOSITE  INTERLAMINAR 
FRACTURE  TOUGHNESS  USING  THE  END-NOTCH-FLEXURE  (ENF) 

OR  MIXED-MODE-FLEXURE  (MMF)  METHOD 


o  DETERMINE  REQUISITE  PROPERITES  AT  DESIRED  CONDITIONS  USING  COMPOSITE 
MICROMECHANICS 


o  RUN  3-D  FINIIt  ELEMENT  ANALYSIS  ON  ENF  (MMF)  FOR  AN  ARBITRARY  LOAD 

«  SCALE  LOAD  TO  MATCH  INTERLAMINAR  SHEAR  STRESS  AT  ELEMENT  NEXT  TO 

CRACK-TIP 


WITH  SCALED  LOAD  EXTEND  CRACK  AND  PLOT  STRAIN  piERGY  RELEASE  RATE 
VS  CRACK  LENGTH  q| 

SELECT  CRITICAL  "G" 

AND  CRITICAL  "a" 


o  METHOD  HAS  VERSATILITY/GENERALITY 


Extended  crack  length,  a 


FIBER  VOLUME  RATIO  EFFECT  ON  INTERLAMINAR 
FRACTURE  TOUGHNESS  PARAMETERS 

END-NOTCH-FLEXURE  (AS/E) 


INTERLAMINAR 

SHEAR 

STRENGTH, 

ksi 


O  INTERLAMINAR  SHEAR  STRENGTH 
A  INTERPLY  LAYER  THICKNESS 


INTERPLY 

LAYER 

THICKNESS, 

in 
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CODSTRAN  RESULTS 


N00[S  446 

OtCKltS  OF  FRtLDOM  8/1 


NODES  446 

DEGREES  Of  FREEDOM  871 


SUCCESSIVE  DAMAGE  EXTENT  AND  DEFECT  GROWTH 
lilSlg  GRAPIIITE/EPOXY  SOLID  LAMINATE 


SUCCESSIVE  DAMAGE  EXTENT  AND  DEFECT  GROWTH 
I+lSIg  GRAPHITE/EPOXY  LAMINATE  WITH  THROUGH-SLIT 


ORIGINAL  MESH :  EUEMENTS  534 

NODES  442 

DECREES  OF  FREEDOM  863 


SUCCESSIVE  DAMAGE  EXTENT  AND  DEEECT  GROWTH 
(±151 3  GRAPHITE/EPOXY  LAMINATE  WITH  THROUGH-HOLE 


FRACTURE  MODES*  OF  [+0J5  6/E  LAMINATES 
(PREDICTED  BY  CODSTRAN) 


PLY  ORIENTATION; 

C+eJs;  e 

IN  DEGREES 

NOTCH  TYPE 

0 

3 

5 

10 

15 

30 

45 

60 

75 

90 

UNNOTCHED- 

LT 

LT 

LT 

LT 

I 

S 

I 

TT 

TT 

TT 

SOLID 

S3 

S3 

S3 

S 

S 

I-^ 

NOTCHED- 

s‘ 

s' 

s' 

S 

S 

I  4 

i-i 

TT 

TT 

TT 

THRU  SLIT 

LT 

LT 

LT 

S 

S 

,4 

4 

1 

4 

NOTCHED- 

S^ 

S 

s 

I 

S 

I 

TT 

TT 

THRU  HOLE 

LT 

LT 

LT 

LT 

S 

TT 

TT 

*  LT  =  LONGITUDINAL  TENSION 
TT  '  TRANSVERSE  TENSION 
S  *  INTRAPLY  SMEAR:  1) 

2) 


INITIAL  FRACTURE  DUE  TO  INTRAPLY  SHEAR  IN  THE  NOTCH  TIP  ZONE 
MINIMAL  INTRAPLY  SHEARING  DURING  FRACTURE 
SOME  INTRAPLY  SHEAR  OCCURRING  NEAR  CONSTRAINTS  (GRIPS) 


4)  OE1.AHINATIONS  OCCUR  IN  NOTCH  TIP  ZONE  PRIOR  TO  ANY  INTRAPLY  DAMAGE 


I  =  INTER PLY  OELAHINATION 
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COMPOSUE  MICROMECHANICS  •  SCALE  LEVELS 


COMPUTER  PROGRAMS  FOR  OOMPOSHES 


0  ICATJ 

0  COBSIRAT^ 

0  CODSTRAN 

0  N.  L.  COBSTRAN 

0  STAEBL 

0  STAT 


INTEGRATED  COMPOSITES  ANALYZER  -  INTERACTIVE 
COMPOSITE  BLADE  STRUCTURAL  ANALYSIS  -  STAND  ALONE 
COMPOSITE  DURABILITY  STRUCTURAL  ANALYSIS  -  STAND  ALONE 
NONLINEAR  COMPOSITE  BLADE  STRUCTURAL  ANALYSIS  -  STAND  ALONE 
STRUCTURAL  TAILORING  OF  ENGINE  BLADES  -  HIGH  TEMPERATURE 
STRUCTURAL  TAILORING  OF  ffiVANCED  TURBOPROPS 


ICAN:  INTEGRATED  COMPOSITES  ANALYZER 


PLY 

COMPOSUE  i 

miCRUMECItANlCS  j 
IHEORY  / 

/ 

/ 

...  / 

lUP-DOlVN 

MAItRIAl.  PROPtHIlES 


P  (0,  T.  M) 


IRACf  U 
OR 

'•DtCOMPOSIIION" 


SIltUCIURAL  lAILORlNG  OF  ENGINE  BLADES  (STAEBL) 


CD'RS-lf.SfiS 
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STRUCTURALLY  TAILORED  SHROUDLESS  BLADES 
WITH  COMPLEX  INTERNAL  STRUCTURE 


- - 

1  PERCENT  CHANGE  FROM  REFERENCE  DES.IG_n - - - 

PARAMETER 

HOLLOW  BLADE  WITH 

SUPERHYBRID 

COMPOSITE  INLAYS 

•  COMPOSITE  _ 

BLADE  CHORD 

-18 

-9 

BLADE  WEIGHT 

-52  , 

-37 

OIRECl  OPERATIONAL 

COST  PLUS  INTEREST 

(DOC  +1) 

.  ENGINE  WEIGHT 

-0.33 

-0.25 

.  ENGINE  COST 

-0.15 

-0.18 

.  MAINTENANCE  COST 

40.03 

40.05 

total' 

-0.N5 

-0.36  _ 

DESIGN  RLOUIREMENTS  SATISFIED: 

.  RESONANCE  MARGINS  {FIRST,  SECOND,  THIRD,  AND  FOURTH  MODES) 
.  FLUTTER-LOG  DECREMENT  (FIRST,  SECOND,  AND  THIRD  MODES) 

.  BIRD  INGESTION  (LOCAL  AND  ROOT  STRESSES) 

.  STEADY-STATE  STRESS  (THROUGHOUT  THE  BLADE) 


OP n nun  blade  design 

EFFECT  OF  IMERnAL  AND  PRESSURE  LOADS 


OPTIMUM  DESIGN  PRESSURE  LOAD  ONLY  PRESSURE  AND  IHERMAL  LOADS 

TFMPERATune  -  ucrmnENT 
PROPERTIES 


WEIGHT  (LBS) 

PERCENT  SPAN 

THICKNESS  (IN) 

CHORD  (IN) 

TIIICK/CHORD 

9.88 

0,  50.  100. 

,97  .08  .10 

3.35  3.66  9.22 

.19  .02  .02 

9.69 

0.  50.  100. 

.98  .08  .08 

3.27  3. GO  9.15 

.19  .02  .02 

CONSTRAINTS 

RESONANCE  MARGINS 

MODE  1 

.05 

.05 

MODE  2 

1.55 

1.97 

MODE  3 

1.71 

1.67 

FLUTTER  CONSTRAINT 

.510 

.520 

ROOT  STRESS 

.782 

.813 

OPIIHIH  BLADE  DESIGN 
EFFECl  OF  IHERMAL  ENVIRONMENl 


OPTIHUri  DESIGN 


CON1RUL  CASE  I 
No.liiERNAL  Loads 


THERMAL  LOADS 
ICHPERAIUBE  iHOEPENOENT 


THERMAL  LOADS 

ICHPERATURE  DEPENDENT 


Weight  (ids) 

PCnCENT  SPAN 

lllICKNESS  (IN) 

Chord  (in) 

Thick/chord 

7.95 

0.  50.  100. 

.29  .13  .08 

1.70  1.80  2.16 

.19  .07  .09 

8,22 

0.  50.  too. 

.29  .15  .08 

1.68  1.65  2.13 

.19  .00  .09 

8,29 

0.  50.  100. 

.29  .15  .08 

1.69  l.BG  2.19 

.19  .08  .09 

CONS TRAIN IS 

Resonance  Margins 

Root  1 

.05 

.06 

.06 

Root  2 

1.27 

1.28 

1.27 

Root  3 

1.G5 

1.67 

1.65 

FivTTER  Constrains 

.988 

.999 

.999 

Root  Stress 

.055 

.062 

.062 
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COMPOSITE  TURBOPROP  BLADE 


COMPOSITE  TURBOPROP  STRESS  ANALYSIS 


BlWl  INKUm  SIRUCIUn 


COMPOSITE  TURBOPROP  VIBRATION  ANALYSIS 


CONSTRAINED  LAYER  DAnRIHC  OF  ADVANCED  PROPFANS 


VIBRATION  MODE 
COMPOSITE  TURBOPROP  BLADE 


ANALYSIS  test 


BENDING 

MEMBRANE 

MEMBRANE 

BENDING 

BENDING 

TRAN.  SHEAR 

1st  BENDING 

154 

158 

147 

156 

1st  EDGE 

338 

339 

334 

325 

2nd  BENDING 

373 

369 

364 

371 

1st  TORSION 

685 

571 

561 

536 

3rd  BENDING 

655 

650 

635 

623 

NATURAL  FREQUENCIES  OF  SR7A 


HIGH  TB'IPERAIURE  COM  POSHES 


TEST  METHODS  AND  CHARACTERIZATION 
0  COMPOSITE  BURNER  LINER 

0  FRS  FOR  SSME  TURBOPUMP  BLADES 


CONaUSIONS 


CURRENT  LEVJIS  RESEARCH  ACTIVITIES  ON  COMPOSITE  MECHANICS /RELATED  AREAS  INCLUDE: 

0  COMPOSITE  MECHANICS,  COMPUTER  PROGRAMS  FOR  arMPOSITES.  HIGH  TEMPERATURE 
COMPOSITES  AND  CaiPOSITE  ENGINE  STRUCTURAL  COMPONENTS 

0  RECENT  PROGRESS  INaUDES; 

0  SIMPLIFIED  DESIGN  ffiOCEDUTE  FOR  GENERAL  LAMINATES 

0  WIDTH,  LOADING  CONDITIONS  AND  ENVIRONMENTAL  EFFECTS  ON  FREE  EDGE  STRESSES 
0  INTERLAMINAR  FRACTURE  IN  COMPOSHES 
0  PROGRESSIVE  FRACTURE  OF  COMPOSITES  WITH/WnHOUT  DEFECTS 
0  CONTINUING  DEVELOPMENT  OF:  ICAN,  CODSTRAN,  N.  L.  COBSTRAfT,  STAEBL,  STAT 
0  INHIATION  OF  RESEARCH  IN  HIGH  TEMPERATURE  COMPOSITES 
0  STRUCTURAL  TAILORING  OF  COMPOSHE  FAN  BLADES,  TURBOPROPS,  SSME  BLADES 
0  THERMOVISCOPLASTIC  STRUCTURAL  AfJALYSIS  OF  TURBINE  BLADES  MADE  FROM 
TUNGSTEN -FIBER  REINFORCED  SUPERALLOYS  -  DEDICATE  ANALYSIS 


APPENDIX  A:  PROGRAM  LISTINGS 


A-1 


AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 
MATERIALS  LABORATORY 

INHOUSE 


ADVANCED  COMPOSITES 

WORK  UNIT  DIRECTIVE  (WUD)  NUMBER  45 

77  April  -  85  April 

WUD  Leader;  James  M.  Whitney 

Materials  Laboratory 

Air  Force  Wright  Aeronautical  Laboratories 
AFWAL/MLBM 

Wright-Patterson  AFB,  OH  45433-6533 
(513)  255-2340  Autovon:  785-2340 

Objective:  The  objective  of  the  current  thrust  under  this  work  is  to  develop  and 

demonstrate  concepts  of  damage  resistance  as  applied  to  fiber  reinforced 
composite  laminates.  Short  term  objectives  (1-2  yrs)  include  the 
following: 

(a)  Development  of  failure  mode  models  with  emphasis  on  delamination  and 
matrix  cracking. 

(b)  Assess  the  role  of  matrix  toughness  in  composite  failure  processes. 

(c)  To  develop  concepts  of  interface/interphase  strengthening. 


CONTRACTS 

IMPROVED,  DAMAGE  RESISTANT  COMPOSITE  MATERIALS 

F33615-84-C-5070 

1  Sep  84  -  1  Feb  88 

Project  Engineer:  James  M.  Whitney 

Materials  Laboratory 

Air  Force  Wright  Aeronautical  Laboratories 
AFWAL/MLBM 

Wright-Patterson  AFB,  OH  45433-6533 
(513)  255-2340  Autovon:  785-2340 

Principal  Investigator:  Ron  Servais 

University  of  Dayton  Research  Institute 
300  College  Park  Avenue 
Dayton,  OH  45469 

Objective:  The  objective  of  this  program  is  to  investigate  from  both  an  experimental 
and  analytical  standpoint  the  potential  of  new  and/or  modifications  of 
existing  polymeric  materials  and  reinforcement  forms  for  use  in  advanced 
composite  materials,  including  processing/mechanical  property  relation¬ 
ships.  Such  materials  are  subsequent  candidates  for  use  in  advanced 
aircraft  and  aerospace  structural  applications. 


3- DIMENSIONAL  RESPONSE  OF  COMPOSITES 
F33615-85-C-5034 
1  Jun  85-1  Dec  87 

Project  Engineer:  Nicholas  J.  Pagano 

Materials  Laboratory 

Air  Force  Wright  Aeronautical  Laboratories 
AFWAL/MLBM 

Wright-Patterson  AFB,  OH  45433-6533 
(513)  255-6762  Autovon:  785-6762 


A-2 


Principal  Investigator:  Sora  R.  Soni 

Adtech  Systems  Research  Inc. 

211  N.  Broad  Street 
Fairborn,  OH  45324 

Objective:  The  objective  of  this  program  is  to  develop  3-dimensional  analytical 

models  capable  of  predicting  the  mechanical  response  of  thick  laminated 
composites  and  residual  stress  failure  of  carbon-carbon  composites. 

INITIAL  IMPACT  DAMAGE  OF  COMPOSITES 

F33615-84-C-5096 

28  Sep  84  -  15  Sep  86 

Project  Engineer:  James  M.  Whitney 

Materials  Laboratory 

Air  Force  Wright  Aeronautical  Laboratories 
AFWAL/MLBM 

Wright-Patterson  AFB,  OH  45433-6533 
(513)  255-2340  Autovon:  785-2340 

Principal  Investigator:  Jonathan  Goering 

Materials  Sciences  Corporation 
Gwynedd  Plaza  II 
Bethlehem  Pike 
Springhouse,  PA  19477 

Objective:  The  objective  of  this  program  is  to  develop  an  analytical  model  for 

predicting  damage  initiation  in  laminated  fiber  reinforced  composites 
subjected  to  low  velocity  impact. 


FUNDAMENTAL  MATRIX  STIFFNESS  FORMULATIONS  FOR  LAMINATE  STRUCTURES 

F33615-83-C-5076 

1  Jun  83  -  31  Mar  86 

Project  Engineer:  Steven  L.  Donaldson 
Materials  Laboratory 

Air  Force  Wright  Aeronautical  Laboratories 
AFWAL/MLBM 

Wright-Patterson  AFB,  OH  45433-6533 
(513)  255-3068  Autovon:  785-3068 

Principal  Investigator:  Henry  T.  Yang  . 

School  of  Aeronautical  §  Astronautical  Engineering 

Purdue  University 

West  Lafayette,  IN  47907 

(317)  494-5117 

Objective:  This  program  will  develop  the  mathematical  formulation  of  the  stiffness 

and  mass  matrices  of  laminated  plates  and  beams,  to  ultimately  obtain 
the  stress  fields,  the  vibrational,  and  the  buckling  response  of 
structural  laminates.  A  through-the-thickness  calculation  shall  be  made 
of  individual  ply  stresses  and  strength  ratios.  The  elements  will  be 
formulated  in  such  a  way  that  they  can  be  simply  implemented  on  micro 
and  minicomputers. 


FAILURE  ANALYSIS  FOR  COMPOSITE  STRUCTURE  MATERIALS 

F33615-84-C-5010 

Jun  84  -  Nov  86 

Project  Engineer:  Frank  Fechek 

Materials  Laboratory 

Air  Force  Wright  Aeronautical  Laboratories 
AFWAL/MLSE 

Wright-Patterson  AFB,  OH  45433-6533 
(513)  255-7483  Autovon:  785-7483 
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Principal  Investigator:  Brian  Smith 

The  Boeing  Company 
P.O.  Box  3707 
Mail  Stop  73-43 
Seattle,  WA  98124 

Objective:  The  objectives  of  this  program  are:  a)  to  verify  the  capability  of  state- 

of-the-art  analysis  techniques  and  procedures  to  produce  useful  data 
toward  the  understanding  o£  the  cause  o£  composite  £ailures,  beginning 
with  the  £ailed  part  and,  b)  to  organize  this  in£ormation  into  a 
compendium  de£ining  a  £ailure  logic  network  which  will  assist  the  failure 
analyst  in  sequentially  selecting  the  appropriate  tests,  techniques,  and 
procedures  to  be  applied  when  conducting  a  failure  analysis  of  a  composite 
structure. 


CURING  PROCESS  OF  COMPOSITE  MATERIALS 

F33615-84-C-5049 

Sep  84  -  1  Oct  87 

Project  Engineer:  Stephen  W.  Tsai 

Materials  Laboratory 

Air  Force  Wright  Aeronautical  Laboratories 
AFWAL/MLBM 

Wright-Patterson  AFB,  OH  45433 
(513)  255-3068  Autovon:  785-3068 

Principal  Investigator:  George  S.  Springer 

Dept  of  Aeronautics  and  Astronautics 
Stanford  University 
Stanford,  CA  94305 
(415)  497-4135 

Objective:  To  extend  the  analytical  modeling  developed  by  the  Principal  Investigator 

to  include  the  curing  thermosetting  and  thermal  plastics  as  the  matrix 
materials.  To  provide  criteria  for  automated  process  controls  and 
optimization. 


AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH 
INHOUSE 


NONE 


GRANTS  AND  CONTRACTS 


DAMAGE  MODELS  FOR  CONTINUOUS  FIBER  COMPOSITES 
84  February  01  -  87  January  31 

Principal  Investigator;  Dr  David  H  Allen 

Department  of  Aerospace  Engineering 
Texas  A&H  University 
College  Station,  TX  77843 
(409)  845-7541 

Program  Manager:  HaJ  David  A  Glasgow 
AFOSR/NA 

Bolling  AFB  DC  20332-6448 
(202)  767-4937 

Objective:  To  develop  a  damage  model  for  predicting  strength  and  stiffness  of  continuous  fiber 

composite  structure  subjected  to  fatigue  loading,  and  to  verify  this  model  with  experimental  results. 


DELAHINATION  AND  TRANSVERSE  FRACTURE  IN  GRAPHITE/EPOXY  COMPOSITES 
84  February  01  -  86  March  31 

Principal  Investigator:  Or  Walter  L  Bradley 

Department  of  Mechanical  Engineering 
Texas  A&M  University 
College  Station,  TX  77843 
(409)  845-1259 

Program  Manager:  MaJ  David  A  Glasgow 
AFOSR/NA 

Bolling  AFB  DC  20332-6448 
(202)  767-4937 

Objective;  To  better  define  the  deformation  and  fracture  physics  of  delamlnatlon  and  transverse 
fracture  In  graphite  epoxy  composites,  and  to  Incorporate  more  realistic  macroscopic  measures  of  the 
fracture  process  Into  linear  and  nonlinear  materials  characterization. 


DYNAMICS  AND  AEROELASTICITY  OF  COMPOSITE  STRUCTURES 
84  May  01  -  86  June  30 

Principal  Investigator:  Dr  John  DugundJI 

Department  of  Aeronautics  &  Astronautics 
Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139 
(617)  253-3758 

Program  Manager;  Dr  Anthony  K  Amos 
AFOSR/NA 

Bolling  AFB  DC  20332-6448 
(202)  767-4937 

Objective;  To  pursue  combined  experimental  and  theoretical  Investigations  of  aeroelastic  tailoring 
effects  on  flutter  and  divergence  of  aircraft  systems. 

ANALYTICAL  AND  EXPERIMENTAL  CHARACTERIZATION  OF  DAMAGE  PROCESSES  IN  COMPOSITE  LAMINATES 
84  September  30  -  85  September  29 


A-5 


Principal  Investigator:  Dr  George  J  Dvorak 

Department  of  Civil  Engineering 
Rensselaer  Polytechnic  Institute 
Troy.  NY  12181 
(518)  266-6943 


Program  Manager:  MaJ  David  A  Glasgow 
AFOSR/NA 

Bolling  AFB  DC  20332-6448 
(202)  767-4937 


Oblective:  To  develop  distributed  damage  analysis  applicable  to  high  matrix  crack  densities,  examine 
damage  propagation  across  and  along  ply  Interfaces,  model  damage  growth 

regions,  and  analyze  stability  and  compressive  strength  of  laminated  plates  containing  distributed 
and/or  concentrated  damage. 


BEHAVIOR  OF  FIBRE  REINFORCED  COMPOSITES  UNDER  DYNAMIC  TENSION 
84  March  15  -  86  March  14 

Principal  Investigators:  Dr  John  Harding 

Dr  C  Ruiz 

Department  of  Engineering  Science 
University  of  Oxford 
Oxford,  0X1  3PJ  England 


Program  Manager:  Or  Anthony  K  Amos 
AFOSR/NA 

Bolling  AFB  DC  20332-6448 
(202)  767-4937 


Objective:  To  characterize  the  mechanical  behavior  and  failure 

Kevlar/epoxy,  and  hybrid  composites  under  tensile  Impact  loading  using  specially  designed  split 
Hopkinson  bar  equipment. 


ANALYSIS  OF  FATIGUE  DAMAGE  AND  FAILURE  IN  COMPOSITE  MATERIALS 
84  September  30  -  87  December  31 

Principal  Investigator:  Or  Zvi  Hashin 

Dept  of  Materials  Science  and  Engineering 
University  of  Pennsylvania 
Philadelphia,  PA  19104 
(215)  898-8337 

Program  Manager:  Ma)  David  A  Glasgow 
AFOSR/NA 

Bolling  AFB  DC  20332-6448 
(202)  767-4937 


Oblective:  To  evaluate  stiffness  change  In  laminates  due  to  distribution  of  ?hf 
Interlaminar  cracks  by  the  use  of  variational  methods,  and  to  determine  the  relationship  between  the 
stiffness  deterioration  and  the  strength  deterioration  of  cracked  laminates. 


RESISTANCE  CURVE  APPROACH  TO  PREDICTING  RESIDUAL  STRENGTH  OF  COMPOSITES 
84  August  01  -  86  Duly  31 

Principal  Investigator:  Dr  H  P  Kan 

Northrop  Corporation 
One  Northrop  Avenue 
Hawthorne,  CA  90250 
(213)  970-2134 

Program  Manager:  Maj  David  A  Glasgow 
AFOSR/NA 

Bolling  AFB  DC  20332-6448 
(202)  767-4937 
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Objective:  To  experimentally  determine  the  Mode  II  delamlnatlon  growth  resistance  of  composite 
laminates  and  to  develop  analytical  techniques  for  application  of  the  R-curve  concept  to  residual 
strength  prediction  of  composite  laminates  with  delamlnatlons . 


ULTRASONIC  NOE  OF  DAMAGE  IN  CONTINUOUS  FIBER  COMPOSITES 
84  February  01  -  87  January  31 

Principal  Investigator:  Dr  Vikram  K  Kinra 

Department  of  Aerospace  Engineering 
Texas  A&H  University 
College  Station,  TX  77843 
(409)  845-1667 

Program  Manager:  Maj  David  A  Glasgow 
AFOSR/NA 

Bolling  AFB  DC  20332-6448 
(202)  767-4937 

Objective:  To  develop,  test,  and  Implement  ultrasonic  nondestructive  evaluation  techniques  to 

characterize  damage  states  produced  In  continuous  fiber  composites  by  monotonic  and  fatigue  loading. 


FRACTURE  AND  LONGEVITY  OF  COMPOSITE  STRUCTURES 
82  January  01  -  85  March  14 

Principal  Investigator:  Or  Paul  A  Lagace 

Department  of  Aeronautics  &  Astronautics 
Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139 
(617)  253-2426 

Program  Manager:  MaJ  David  A  Glasgow 
AFOSR/NA 

Bolling  AFB  DC  20332-6448 
(202)  767-4937 

Objective:  To  develop  theoretical  and  seml-empirical  fracture  laws  and  failure  criteria  and  to 

correlate  them  with  extensive  experimental  data  generated  In  the  program. 


NONLINEAR  DYNAMIC  RESPONSE  OF  COMPOSITE  ROTOR  BLADES 
82  September  01  -  86  November  30 

Principal  Investigators:  Dr  Ozden  Ochoa 

Department  of  Mechanical  Engineering 
Texas  A&M  University 
College  Station,  TX  77843 
(409)  845-2022 


Or  John  J  Engblom 

Department  of  Mechanical  Engineering 
Texas  A&M  University 
College  Station,  TX  77843 
(409)  845-2813 

Program  Manager:  Dr  Anthony  K  Amos 
AFOSR/NA 

Bolling  AFB  DC  20332 
(202)  767-4937 

Objective:  To  develop  nonlinear  finite  element  models  suitable  for  predicting  the  structural  dynamic 
response  and  resulting  damage  of  composite  rotor  blades  under  Impact  and  other  transient  excitations. 


INTERLAMINAR  FRACTURE  TOUGHNESS  IN  RESIN  MATRIX  COMPOSITES 
83  January  01  -  87  April  14 
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Principal  Investigator:  Dr  Lawrence  W  Rehfleld 

School  of  Aerospace  Engineering 
Georgia  Institute  of  Technology 
Atlanta,  6A  30332 
(404)  894-3067 


Program  Manager:  MaJ  Oavid  A  Glasgow 
AFOSR/NA 

Bolling  AFB  DC  20332-6448 
(202)  767-4937 


Objective:  To  develop  a  Mode  II  Interlaminar  fracture  coupon  and  test  that  can  be  used  tn  both 

tension  and  compression  testing,  can  be  analyzed  ^  ^ 

Interpreted  and  provides  an  experimental  means  for  Isolating  Mode  II  contributions  to  fracture. 


DAMAGE  MODELS  FOR  DELAMINATION  AND  TRANSVERSE  FRACTURE  IN  FIBROUS  COMPOSITES 
84  February  15  -  86  February  14 

Principal  Investigator:  Or  Richard  A  Schapery 

Department  of  Civil  Engineering 
Texas  A&M  University 
College  Station,  TX  77843 
(409)  845-7512 


Program  Manager:  Maj  David  A  Glasgow 
AFOSR/NA 

Bolling  AFB  DC  20332-6448 
(202)  767-4937 


Objective:  To  develop  and  verify  mathematical  models  of  delamlnatlon  and  transverse  fracture  which 
Sount  for  local  (crack  tip)  and  global  damage  distributions,  separating  the  lay-up  dependent 
fracture  energy  associated  with  microcracking  from  the  Intrinsic  fracture  energy  of  the  separation 
.process  which  occur  at  the  tip  of  an  advancing  delamlnatlon  crack. 


EFFECT  OF  LOCAL  MATERIAL  IMPERFECTIONS  ON  BUCKLING  OF  COMPOSITE  STRUCTURAL  ELEMENTS 
83  June  30  -  85  August  31 

Principal  Investigator:  Dr  George  J  Simitses  ,  „  .  . 

Dept  of  Engineering  Science  and  Mechanics 
Georgia  Institute  of  Technology 
Atlanta,  GA  30332 
(404)  894-2770 

Program  Manager:  Dr  Anthony  K  Amos 
AFOSR/NA 

Bolling  AFB  DC  20332-6448 
(202)  767-4937 


Objective:  To  Investigate  the  effects  of 
the  buckling  characteristics  of  composite 


localized  material,  geometric,  and  process  Imperfections  on 
structural  elements,  and  to  Incorporate  them  In  analytical 


prediction  methods. 

COMPREHENSIVE  STUDY  ON  DAMAGE  TOLERANCE 
84  September  30  -  85  September  29 


PROPERTIES  OF  NOTCHED  COMPOSITE  LAMINATES 


Principal  Investigator:  Or  Albert  S  0  Wang 

Dept  of  Mechanical  Engineering  and  Mechanics 
Orexel  University 
Philadelphia,  PA  19104 
(215)  895-2297 


Program  Manager:  Maj  David  A  Glasgow 
AFOSR/NA 

Bolling  AFB  DC  20332-6448 
(202)  767-4937 

Objective:  To  conduct  a  comprehensive  analysis  of  the  stress  fields  In  notched  laminates  so  as  to 
develop  a  fundamental  understanding  of  the  damage  mechanisms  near  the  notch  region. 
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RESIDUAL  STRESS  INDUCED  DAMAGE  IN  COMPOSITE  MATERIALS 
84  February  01  -  85  January  31 

Principal  Investigator;  Dr  Y  Weltsman 

Department  of  Civil  Engineering 
Texas  A&M  University 
College  Station,  TX  77843 
(409)  845-7512 

Program  Manager:  MaJ  David  A  Glasgow 
AFOSR/NA 

Bolling  AFB  DC  20332-6448 
(202)  767-4937 

Objective:  To  develop  and  verify  methods  for  predicting  damage  formation,  growth,  and  arrest  due  to 
residual  stresses  In  fiber-reinforced,  resin  matrix  composites. 
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OFFICE  OF  NAVAL  RESEARCH 
ARLINGTON,  VA  22217 

CONTRACTS 


FLAW  GROWTH  AND  FRACTURE  OF  COMPOSITE  MATERIALS 
AND  ADHESIVE  JOINTS 
N00014-79-C-0579 
July  83  -  Nov  87 

Project  Engineer:  Dr.  Yapa  Rajapakse 

OFFICE  OF  NAVAL  RESEARCH 
Mechanics  Division,  Code  432S 
Arlington,  VA  22217 
(202)  696-4306  Autovon  226-4306 


Principal  Investigator:  Prof.  S.  S.  Wang 

University  of  Illinois 

Department  of  Theoretical  and  Applied  Mechanics 
Urbana,  Illinois  61801 
(217)  333-1835 


Objective: 


Analytical  and  numerical  studies  will  be  conducted  of  flaw  growth  and 
Fracture  in  Fiber  Composite  Laminates  and  adhesively  bonded  structural 
joints  under  static  and  dynamic  loading  conditions . 


DAMAGE  ACCUMULATION  AND  RESIDUAL  PROPERTIES  OF  COMPOSITES 

N00014-82-K-0572 

July  82  -  March  86 

Project  Engineer:  Dr.  Yapa  Rajapakse 

OFFICE  OF  NAVAL  RESEARCH 
Mechanics  Division,  Code  432S 
Arlington,  VA  22217 
(202)  696-4306  Autovon  226-4306 


Principal  Investigator:  Prof.  I.  M.  Daniel 

Illinois  Institute  of  Technology 
Department  of  Mechanical  Engineering 
Chicago,  Illinois  60616 
(312)  567-3186 


Objective; 


Investigate  damage  mechanisms  and  damage  accumulation  in  graphite /epoxy 
laminates  for  the  development  of  models  for  predicting  residual 
stiffness,  residual  strength,  and  residual  life. 


INVESTIGATIONS  OF  ENVIRONMENTAL  EFFECTS  AND  ENVIRONMENTAL 

DAMAGE  IN  COMPOSITES 

N00014-82-K-0562 

October  84  -  September  87 


Project  Engineer:  Dr.  Yapa  Rajapakse 

OFFICE  OF  NAVAL  RESEARCH 
Mechanics  Division,  Code  423S 
Arlington,  VA  22217 
(202)  696-4306  Autovon  226-4306 


Principal  Investigator:  Prof.  Y.  Weitsman 

Texas  A&M  University 
Department  of  Civil  Engineering 
College  Station,  Texas  77843 
(713)  845-7512 


Objective : 


Research  will  be  conducted  to  study  the  effects  of  stress  and  moisture  on 
the  mechanical  response  of  graphite /epoxy  composites.  Special  attention 
V7ill  be  given  to  environmental  induced  damage  growth  and  its  effect  on 
compressive  and  shear  response. 
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INVESTIGATIONS  OF  IMPACT  DAMAGE  IN  COMPOSITES 

N00014-84-K-0460 

June  84  -  Jan  88 

Project  Engineer:  Dr,  Yapa  Rajapakse 

OFFICE  OF  NAVAL  RESEARCH 
Mechanics  Division,  Code  423S 
Arlington,  VA  22217 
(202)  696*“4306  Autovon  226-4306 

Principal  Investigator:  Prof.  J.  Awerbuch 

Drexel  University 

Department  of  Mechanical  Engineering  and  Mechanics 
Philadelphia,  PA  19104 
(215)  895-2291 

Objective:  Investigations  of  damage  in  graphite /epoxy  laminates  due  to  normal  and 

oblique  impact  will  be  carried  out  using  a  variety  of  experimental 
techniques.  The  use  of  acoustic  emission  for  damage  assessment  will  be 
explored  fully. 

SUPPRESSION  OF  DELAMINATION  IN  COMPOSITE  LAMINATES  SUBJECTED 
TO  IMPACT  LOADING 
N00014-84-K-0554 
July  84  -  June  86 

Project  Engineer:  Dr.  Yapa  Rajapakse 

OFFICE  OF  NAVAL  RESEARCH 
Mechanics  Division,  Code  423S 
Arlington,  VA  22217 
(202)  696-4306  Autovon  226-4306 


Principal  Investigator:  Prof.  C.  T.  Sun 

Purdue  University 
West  Lafayette,  IN  47907 
(317)  494-5130 

Obiective:  Research  will  be  performed  to  investigate  and  establish  quantative  models 

for  delamination  growth  in  composite  laminates  specif ically  designed  to 
suppress  de lamination  by  the  use  of  3-D  stitching  reinforcement  and  by 
the  introduction  of  soft  adhesive  layers. 

CONSTRUCTION  OF  NON-LINEAR  MODEL  FOR  BINARY  METAL  MATRIX  COMPOSITES 

N00014-84-K-0468 

July  84  -  June  86 

Project  Engineer:  Dr.  A.  S.  Kushner 

OFFICE  OF  NAVAL  RESEARCH 
Mechanics  Division,  Code  423S 
Arlington,  VA  22213 
(202)  696-4306  Autovon  226-4306 


Principal  Investigator:  Prof.  H.  Murakami 

University  of  California,  San  Diego 
La  Jolla,  CA  92093 
(619)  452-3821 


Obiective:  Non-linear  theory  for  metal  matrix  composites  will  be  developed,  based  on 

^  variational  principles  and  multi-variable  asmptotic  expansion  techniques. 

The  theory  will  account  for  the  effect  of  fiber  breakage,  fiber-matrix 
debonding  and  slip,  matrix  plasticity  and  delamination. 


METAL  MATRIX  COMPOSITE  INTERFACES 

N00014-84-K-0495 

September  84  -  August  87 

Project  Engineer:  Dr.  A.  S.  Kushner 

OFFICE  OF  NAVAL  RESEARCH 
Mechanics  Division,  Code  423S 
Arlington,  VA  22213 
(202)  692-4306  Autovon  226-4306 
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Principal  Investigator:  Prof.  A.  S.  Argon 

MIT 

Department  of  Mechanical  Engineering 
Cambridge,  MA  02139 
(617)  253-2217 

Objective:  Research  will  be  conducted  to  develop  the  micro-mechanical  model  of  the 

interface  in  metal  matrix  composites  which  have  the  features  of  predict¬ 
ability  for  the  purpose  of  optimizing  existing  fiber-matrix  systems. 

FRACTURE  OF  FIBROUS  COMPOSITES  AND  LAMINATES 

N00014-85-K-0247 

March  85  -  Apr  87 

Project  Engineer:  Dr.  Yapa  Rajapakse 

OFFICE  OF  NAVAL  RESEARCH 
Mechanics  Division,  Code  423S 
Arlington,  VA  22217 
(202)  696-4306  Autovon  226-4306 

Principal  Investigator:  Prof.  G.  J.  Dvorak 

Rensselaer  Polytechnic  Institute 
Department  of  Civil  Engineering 
Troy,  NY  12181 
(518)  266-6363 

Objectives  Research  will  be  conducted  to  investigate  the  fracture  behavior  of 
notched  metal  matrix  composites,  accounting  for  plasticity  effects. 

MECHANICAL  PROPERTIES  OF  COMPOSITES  AT  ELEVATED 

TEMPERATURES 

N00014-85-K-0480 

July  85  -  June  88 

Project  Engineer:  Dr.  Yapa  Rajapakse 

OFFICE  OF  NAVAL  RESEARCH 
Mechanics  Division,  Code  423S 
Arlington,  VA  22217 
(202)  696-4306  Autovon  226-4306 

Principal  Investigator:  Prof.  G.  S.  Springer 

Stanford  University 

Dept,  of  Aeronautics  and  Astronautics 
Stanford,  CA  94305 
(415)  497-4135 

Objective:  Mechanics-based  models  for  changes  in  mechanical  properties  and  failure 

characteristics  of  organic— matrix  composites  will  be  established. 

PRECISION  ULTRASONIC  MEASUREMENTS  IN  COMPOSITE 

MATERIALS 

N00014-85-K-0595 

July  85  -  June  87 

Project  Engineer:  Dr,  Yapa  Rajapakse 

OFFICE  OF  NAVAL  RESEARCH 
Mechanics  Division,  Code  423S 
Arlington,  VA  22217 
(202)  696-4306  Autovon  226-4306 

Principal  Investigator:  Prof.  W,  Sachse 

Cornell  University 

Dept,  of  Theoretical  and  Applied  Mechanics 
Ithaca,  NY  14853 

Objective:  Research  will  be  conducted  to  establish  the  capability  to  distinguish 

between  the  different  failure  modes  in  composite  by  the  analysis  of 
acoustic  emission  signals. 
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NAVAL  AIR  SYSTEMS  COMMAND 
WASHINGTON,  D.C.  20361 

INHOUSE 

DEVELOPMENT  OF  CERTIFICATION  METHODOLOGY  FOR 
COMPOSITE  MATERIALS 
October  84  ~  October  87 

Principal  Investigator:  Prof.  E.  Wu 

Naval  Post  Graduate  School 
Department  of  Aeronautics,  Code  67WT 
Montery,  CA  93943 
(408)646-3459  Autovon  878-3459 

Objective:  Develop  probabilistic-based  certification  methods  through  experimental  and 

probabilistic  modeling  to  insure  the  strength  and  life  of  critical 
composite  structure. 


CONTRACTS 

FATIGUE  LIFE  AND  RESIDUAL  STRENGTH  OF  COMPOSITE  STRUCTURES 
September  83  -  September  85 

Project  Engineer:  Dr.  D.  R.  Mulville 

Naval  Air  Systems  Command 
Washington,  D.C.  20361 
(202)  692-7443  Autovon  222-7443 

Principal  Investigators:  Dr.  J.  Yang  and 

Dr.  D.  Jones 

The  George  Washington  University 
Washington,  D.C.  20052 
(202)  676-6929 

Obiective:  Develop  statistical  models  to  describe  fatigue  life  and  residual  strength 

of  composite  structures  including  bolted  and  bonded  composite  joints. 

STRENGTH  AND  IMPACT  RESISTANCE  OF  ADVANCED  COMPOSITE 
SANDWICH  STRUCTURES 
N00019-85-C-0090 
December  84  -  December  87 

Project  Engineer:  Dr.  D.  R.  Mulville 

Naval  Air  Systems  Command 
Washington,  D.C.  20361 
(202)  692-7443 

Principal  Investigators:  Prof.  E.  A.  Witmer  and 

Prof.  P.  A.  Lagace 

Massachusetts  Institute  of  Technology 
Department  of  Aeronautics  and  Astronautics 
Cambridge,  MA  02139 
(617)  253-3628 

Objective:  Further  develop  the  technology  for  advanced  composite  structures  with 

sandwich  or  honeycomb  cores. 

NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 

INHOUSE 

STRUCTURAL  RESPONSE  OF  DAMAGED  COMPOSITES 
October  79  -  September  86 

Project  Engineer:  Dr.  P.  W.  Mast 

Naval  Research  Laboratory 
Washington,  D.C.  20375 
(202)767-2165  Autovon  297-2165 

Objective:  Develop  a  capability  for  predicting  the  structural  response  of  composite 

structures  containing  a  defect  or  damage. 
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NAVAI^  air  PEVELOPWENg?  CENTER 
AIRCRAFT  AND  CREW  SYSTEMS  TECHNOLOGY  DIRECTORATE 
WARMINSTER,  PA  18974 

INHOUSE 


IMPROVED  MATRIX  DOMINATE  PROPERTIES 
October  82  -  September  86 

Project  Engineer:  Ramon  Garcia 

Naval  Air  Development  Center 
ACSTD/6043 

Warminster,  PA  18974 
(215)  441-1321  Autovon  441-1321 

Objective:  Investigate  methods  to  improve  composite  performance  by  reinforcing  the 

matrix  with  silicon  carbide  wiskers. 

HYBRID  COMPOSITE  FRACTURE  CHARACTERIZATION 
September  85  -  October  87 

Project  Engineer:  Lee  W.  Gause 

Naval  Air  Development  Center 
ACSTD/6043 

Warminster,  PA  18974 
(215)  441-1330  Autovon  441-1330 

Objective:  Characterize  the  strength,  mechanical  properties,  and  damage  tolerance 

of  woven  and  hybrid  composite  structures. 


METAL  MATRIX  CRACK  INITIATION/PROPAGATION 
September  85  -  October  87 


Project  Engineer:  Dr.  H.  C.  Tsai 

Naval  Air  Development  Center 
ACSTD/6043 

Warminster,  PA  18974 
(215)  441-2871  Autovon  441-2871 


Obiective:  Characterize  the  crack  initiation/propagation  mechanics  of  silicon  carbide/ 

Objective,  cnarac^^  composites  as  applied  to  landing  gear  and  arrester 

hooks  in  the  Naval  shipboard  environment. 
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DESIGN  OF  HIGHLY  LOADED  COMPOSITE  JOINTS  AND 
ATTACHMENTS  FOR  TAIL  STRUCTURES 
N62269-82-C-0239 
February  82  -  January  86 

Project  Engineer:  Ramon  Garcia 

Naval  Air  Development  Center 
ACSTD/6043 

Warminster,  PA  18974 
(215)  441-1321  Autovon  441-1321 


Principal  Investigator:  S.  W.  Aver ill 

Northrop  Corporation 
Aircraft  Group 
Hawthorne,  CA  90250 
(213)  970-3442 


Obiective:  To  develop  composite  designs  which  will  permit  the  use  of  metal  to 

^  composite  bolted  root  attachments  in  aircraft  tail  structures  as  an 

alternative  to  high-load  transfer  adhesive  bonded  titani^  step  joints. 

To  improve  damage  tolerance,  survivability  and  repairability  over  current 
composite  designs.  Structural  efficiency,  manufacturing  feasibility  an 
quality  assurance  requirements  will  be  determined. 
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DESIGN  OF  HIGHLY  LOADED  COMPOSITE  JOINTS 
AND  ATTACHMENTS  FOR  WING  STRUCTURES 
N62269-82-C-0238 
February  82  -  December  85 

Project  Engineer:  Ramon  Garcia 

Naval  Air  Development  Center 
ACSTD/6043 

Warminster,  PA  18974 
(215)441-1321  Autovon  441-1321 

Principal  Investigator:  M.  J.  Ogonowski 

McDonnell  Aircraft  Co. 

P.O.  Box  516 

St.  Louis,  MO  63166 

(314)  233-8630 

Objective:  To  develop  composite  designs  which  will  permit  the  use  of  metal  to 

composite  bolted  root  attachments  in  aircraft  wing  structures  as  an 
alternative  to  high-load  transfer  adhesive  bonded  titanium  step  joints. 
Strain  concentration  around  fastener  holes,  fatigue  and  environmental 
affects,  damage  tolerance  and  repairability  for  each  concept  will  be 
determined. 

IMPACT  RESPONSE  OF  COMPOSITE  STRUCTURES 

N62269-85-M-3131 

March  85  -  September  85 

Project  Engineer:  Lee  W.  Cause 

Naval  Air  Development  Center 
ACSTD/6043 

Warminster,  PA  18974 
(215)  441-1330  Autovon  441-1330 

Principal  Investigator:  Prof.  P.  V.  McLaughlin 

Villanova  University 

Department  of  Mechanical  Engineering 
Villanova,  PA  19085 
(215)  645-4991 

Objective:  Analytically  describe  the  effects  of  viscoelasticity,  contact  deformation; 

and  shear  deformation  to  the  impact  response  of  composite  material  and 
correlate  these  results  to  observed  experimental  behavior. 

INFLUENCE  OF  LOAD  FACTORS  AND  TEST  METHODS  ON 

IN-SERVICE  RESPONSE  OF  COMPOSITE  MATERIALS  AND 

STRUCTURES 

N62269-85-C-0234 

June  85  -  June  88 

Project  Engineer:  Lee  W.  Cause 

Naval  Air  Development  Center 
ACSTD/6043 

Warminster,  PA  18974 
(215)  441-1330  Autovon  441-1330 

Principal  Investigator:  Prof.  K,  L.  Reif snider 

Virginia  Polytechnic  and  State  University 
Depart,  of  Engineering  Science  &  Mechanics 
Blacksburg,  VA  24061 
(703)  961-5316 

Objective:  Develop  an  understanding  of  the  relationship  between  composite  laminate 

response  to  high  load  levels  for  short  time  periods  and  response  to  low 
load  levels  for  long  time  periods.  Develop  an  understanding  of  the 
relationship  between  test  methods  and  laminate  response.  Establish  the 
manner  in  which  these  relationships  are  associated  with  strength  and  life 
Formulate  a  mechanistic  model  which  can  be  used  to  anticipate  long-term 
behavior . 
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LAMINATED  COMPOSITE  MIXED-MODE  PLANE  CRACK  GROWTH  CRITERIA 

N62269-85-C-0246 

June  85  -  June  86 

Project  Engineer:  Lee  W.  Cause 

Naval  Air  Development  Center 
ACSTD/6043 

Warminster,  PA  18974 
(215)  441-1330  Autovon  441-1330 

Principal  Investigator:  Prof.  A.  S.  D.  Wang 

Drexel  University 

Department  of  Mechanical  Engineering  &  Mechanics 
Philadelphia,  PA  19104 
(215)  895-2297 

Objective:  Define  the  critical  conditions  necessary  to  propagate  a  mixed-mode  plane 

crack  in  composite  laminates  for  use  in  the  automated  crack -growth 
.simulation  program. 
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FORMULATION  OF  LAMINATED  BEAM  AND  PLATE 
FINITE  ELEMENTS  FOR  A  MICROCOMPUTER 
A.  T.  Chen  and  T.  Y.  Yang 
School  of  Aeronautics  and  Astronautics 
Purdue  University 
West  Lafayette,  Indiana  47906 


ABSTRACT 

The  purpose  of  this  paper  is  to  develop  simple  yet  efficient 
formulation  for  laminated  composite  finite  elements  and  also  to 
develop  highly  efficient  numerical  algorithms  for  structural 
analysis  and  design  using  stand  alone  desktop  microcomputers. 

In  the  finite  element  formulations,  an  8  degree  of  freedom 
beam  element  and  an  18  d.o.f  triangular  plate  element  in  bending 
with  anisotropic  symmetrically  laminated  composite  materials  are 
formulated.  In  the  development  of  numerical  procedures,  emphasis 
has  been  placed  upon  minimization  and  condensation  of  memory 
storage  and  efficiency  of  computation  so  that  the  present 
development  is  suitable  for  implementation  on  desktop 
microcomputers.  The  present  development  is  not  only  limited  to 
linear  static  stress  analysis,  it  also  includes  eigenvalue 
analysis  for  free  vibration  and  buckling  problems.  For  the 
special  case  of  lumped  mass  matrices,  a  special  condensation 
technique  making  use  of  zero  terms  along  the  diagonal  of  the  mass 
matrix  is  used. 

To  demonstrate  and  evaluate  the  present  development, 
numerical  analyses  on  the  static,  free  vibration,  and  buckling 
analysis  of  anisotropic  symmetrically  laminated  beam  and  plate 
problems  have  been  analyzed  using  a  microcomputer. 


FIG.  6.  TECHNIdOE  TO  HAKE  EIGENVfiLOE  HATEIK  SyHHETRIC 


Table  1.  liar i oils, ScheHes  for  MeHorn)  Storage  for  [Kl 
for  a  SiMplg  Supported  Square  Isotropic 
Plate  Onder  UniforH  Load 

Kesh  Pattern 


Plate  EleKent  (18  d.o.f,  triangular) 


Mesh  for 
a  Quadrant 

Degrees  of 
FreedoH 

Storage  Modes 
Full  Matrix  SynKetric 

Band 

Sbgline 

3x3 

60 

3600 

1830 

1395 

1001 

4x4 

104 

10816 

5460 

3569 

2?3? 

5x5 

160 

25600 

12880 

6664 

5257 

Table  2,  CoHparison  of  Centerline  Deflections 
of  Two  Modellings  for  a  llniforslg 
Loaded  Cantilever  BeaH 


40  18  d.o.f,  Plate  EleKents 
Five  12  d.o.f,  Bean  Elements 


TaUe  3. 


CoHParison  of  Centerline  Twist 
of  two  Hodellings  for  a  Dnifornly 
Loaded  Cantilever  BeaH 


40  18  d.o.f.  Plate  Elesents 


^  1 1 1 1 '  I  M-q 

Five  12  d.o.f.  BeaH  Elenents 


Fiber 

Angle 

BeaH  Twist 

Tip  x-0.6t  x=0i6L/x=t 

U)  (2)  (2)/(l) 

Plate  Twist 

Tip(x=L)  x=0.6L  x:0.6I/x:i 
(3)  (4)  (3)/(4) 

0 

30 

45 

60 

90 

0.0  0.0 

-3.126  -2.885  0.9231 

-3.047  -2.813  0.9232 

-2.153  -1.987  0.9231 

0.0  0.0 

0.0  0.0 

-2.687  -2.483  0.9240 

-2.787  -2.597  0.9317 

-2.046  -1.912  0.9347 

0.0  0.0 

TaLle  4.  Natural  Frequencies  (Hz)  for  a  Thin  Anisotropic 
6  layer  10  /01s  Cantilever.  Plate 
2 


Layup 

Sequence 


Limped  Kass 
(1)  (2) 

8  el  32  el  8el  32  el 
50  dof  147  dof  8  dof  24  dof 


Consistent  Hass 
50  dof  147  dof 


REF  [11 
Exp  365 
dof 


10.7 

10.9 

10.7 

11.0 

11.0 

[Oj  /981s 

21.0 

31.0 

22.9 

32.2 

39.5 

56.9 

66.4 

59.2 

67.0 

69.2 

8.6 

8.8 

8.8 

8.9 

0.9 

[15 2  /81s 

24.7 

35.2 

26.7 

36.7 

42.9 

47.4 

56.8 

58.3 

57.6 

62.9 

6.1 

6.2 

6.2 

6.3 

6.3 

tSOj  /01s 

27.2 

34.5 

28.9 

35.6 

37.3 

35.4 

45.9 

36.9 

47.1 

58.0 

4.7 

4.8 

4.8 

4.8 

4.9 

[452 /81s 

25.3 

28.8 

27.4 

29.4 

30.1 

27.6 

39.0 

29.0 

40.4 

50.8 

4.1 

4.1 

4.1 

4.2 

4.2 

[682  /01s 

22.7 

25.1 

23.5 

25.5 

26.1 

22.9 

32;8 

24.9 

24.1 

42.8 

3.8 

3.8 

3.8 

3.9 

3.9 

[752  /01s 

19.8 

23.5 

20.9 

23.7 

24.3 

21.5 

28.9 

22.9 

30.1 

37.3 

3.7 

3.8 

3.7 

3.8 

3.8 

[902 /01s 

18.7 

23.0 

19.9 

23.2 

23.8 

21.1 

27.6 

22.3 

28.7 

35.2 

8.9 


6.3 


bf  I  ■  w 

4.9 


4.2 


1  « 
3.9 


WV  I  w 

3.8 


11.2  11.1 

42.4  39.5 

70.5  S9.5 
9.4  8.9 

45.8  42.9 

66.2  (2.7 
(.(  (.3 
40.0  37.3 

59.1  56.9 

4.8  4.9 

29.8  30.1 

51.3  49.4 

4.3  4.2 

27.1  26.1 
47.7  41.7 

3.8  3.9 

25.1  24.3 

38.9  36.7 
3.7  3.8 

24.3  23.9 

38.2  35.1 


Table  5.  CoRparison  of  Hatural  Frequencies  (Hz)  of 
Two  Howell ings  for  a  Cantilever  Beaw 


Layup 

BeaH  EleHent 

Plate  Elenent 

Sequence 

Consistent  Mass 

,  LuHped  Hass 

Consistent  Mass 

4  el 

10  el 

20  el 

40  el 

29  01 

40  el 

17  dof 

52  dof 

20  dof 

30  dof 

122  dof 

183  dof 

0.01026 

iPBI 

0.01022 

0.01022 

0.01026 

0.01026 

[0V 

ms 

0.06436 

0.04751 

0.06301 

0.06432 

0.06432 

0.07716 

0.07716 

0.06322 

0.06755 

0.08235 

0.08234 

0.00525 

0.00525 

0.00549 

0.00549 

0.00551 

0.00551 

[30. 

/Bis 

0.03278 

0.03274 

0.03374 

0.03390 

0.03440 

0.03439 

0.09172 

0.09092 

0.06970 

0.09393 

0.09651 

0.09646 

0.00424 

0.00424 

0.00435 

0.00435 

0.00436 

0.00436 

[45. 

/01s 

0.02655 

0.02651 

0.02688 

0.02690 

0.02730 

0.02730 

0.07456 

0.07404 

0.06260 

0.07455 

0.07645 

0.07645 

0.00378 

0.00378 

0.00382 

0.00382 

0.00383 

0.00383 

[60, 

/01s 

0.02374 

0.02371 

0.02363 

0.02364 

0.02399 

0.02400 

0.06685 

0.06634 

0.05351 

0.06548 

0.06715 

0.06717 

[90, 

0.00353 

0.00353 

0.00352 

0.00352 

0.00354 

0.00354 

/01s 

0.02218 

0.02215 

0.02180 

0.02181 

0.02216 

0.02216 

b 

0.06250 

0.06204 

0.04560 

0.06023 

0.06204 

0.06204 

Table  6.  Critical  Bucbling  Loads  for  a  SiHply 
Supported  Anisotropic  Square. Plate 
llnder  UniforH  CoHpressive  Axial  Load 


■ 

n 

Exp.  [-31 

Exact  [^1 

WM 

318.86 

285 

271 

318.91 

387.90 

425 

399 

“ 

ifl 

354.67 

406 

364 

mm 

355.40 

381 

433 

90 

203.88 

210 

251 

203.75 
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STRENGTH,  DEFLECTIONS  AND  IMPACT  DAMAGE  JOINT  NAVY /FA A  EFFORT 

IN  ADVANCED  COMPOSITE  STRUCTURES 

NAVAIR,  :  STRENGTH  AND  DEFLECTIONS  OF  ADVANCED 
NADC  COMPOSITE  SANDWICH  STRUCTURES 

PAUL  A.  LAGACE 
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20  seconds  CPU  time  for  one  complete  loading 


BUCKLING  JIG 


EXPERIMENTAL  APPARATUS:  LOAD  VS.  CENTER  DEFLECTION 

9-TRANSDUCER  RACK 


Load  =  62  kH 

EA01  LINE  REPRESENTS  0.25  MM  DEFLECTION 


OVERVIEW :  IMPACT  RESPONSE 


Distributed  shear 


LOCAL  CONTACT  PROBLEM  SOLUTION  OVERVIEW  LOAD  VS.  APPROACH 
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Normalized  radius,  r/R. 


GLOBAL  PLATE  BENDING  SOLUTION  OVERVIEW 


o 

</) 

o 

z 

M 

Cl 

z 

UJ 

GQ 

UJ 

< 

-J 

Q. 


CO 

O 

J 

CD 


CO 

3 

I- 

< 

h- 

CO 


“0 

> 

4) 

-L 

— 

> 

T) 

L 

L 

0 

4) 

0 

0 

C 

3 

0) 

0. 

c 

E 

0 

W 

0 

3 

13 

0 

c 

9 

0 

0 

£ 

> 

0 

V) 

V 

L 

X 

L 

U 

— 

a 

CO 

3 

3 

(/) 

_c 

a 

c 

« 

£ 

0 

L 

— 

0 

— 

3 

c 

o 

13 

®  “0 

c  0> 
®  0) 
£  3 
4)  ^ 
—  ^ 
a-0 
E 

“  > 

C  0) 
0 

-  4) 

3  *0. 
—  £ 
0  0 
CO  u 


OD 


Cl 

£ 


L. 

0 

£ 

3 

9 


0 15 

L. 
0 
a 

L  £ 
-  0 
u 

“D 

C  4J 
0  _Q 

tfl  — 

C  — 
D  - 
L  3 

0^  V) 
i.  -*-* 
0  — 
3  D 
W 

Lk  4) 
0 


»-  c 

CL  O 
s:  L 
■D  h- 

co 

CO  I 

< 


Wi  w  w 

-  c  L  a 

r  -  rt  c 


O  U_ 


u  u. 

0  E  - 


V)  V) 

CO  w  c 
0  ^  _ 

3  Z  Q  CO  h-! 
CL 


a- 
a  “ 
3  c 


ID 

o 


o 

S3 


4) 

3  C 
“0  O 


0  <33 

CLU- 

<i3 

C  "0 
O 

L 

0)  0) 
u 

L  C 
0  0) 
LU  O 


4) 

0 

o. 

<V 

V 

3 

H“ 

” 

(A 

X 

UJ 

3 

L 

4) 

<o 

X 

0 

N 

-— 

z 

Tj' 

w 

W 

> 

HD 

4> 

\ 

(A  L 

CD 

C 

E 

3  0 

0) 

0  X 

H- 

U) 

CD 

O 

-*  C 

L  3 

O 

_ 

4J 

0  0 

< 

0 

“D 

>  S3 

a. 

-0 

4) 

u- 

0 

4)  4) 

L 

U 

—  03 

H 

4) 

> 

X  X 

> 

C  4) 

Li. 

4) 

£ 

4) 

"D 

0 

0 

X  W 

O 

0 

0 

3 

0  0 

<o 

LiJ 

•T? 

E 

C 

4) 

E 

ve 

L 

X  0 

M 

E 

4) 

X 

4)  > 

L 

4) 

3 

X 

0  X 

W  u 

L 

3  — * 

.J 

CSJ 

C  4) 

0 

L  - 

M 

4)  — 

u 

^  3 

QQ 

> 

£  — 

4) 

W  « 

< 

c 

0 

0 

>-u 

L 

C  W  C 
0  C  0 

Q- 

— '  (/) 

E 

U  4)  “ 

< 

X 

-  C 

3 

£ 

(_) 

(J 

u  - 

£ 

4)  - 

c 

0 

..  — 

£  0  X 

3 

—  “D 

X  — 

0  03  C 

0 

i>  c 

a  0 

L  a  0 

—1 

>  0 

r  X 

Lu  V3  U 

• 

• 

• 

• 

ONLlf  steel  balls  used 


SPECIMEN  CONFIGURATION  FOR 

PHASE  n  TEST  PROGRAM  FOR  IMPACT  RESPONSE  IMPLANTED  DELAMINATION  TESTS 


as 

-5 

cn  <;> 
u 

s  > 

25 


•  0  O 

ill 


oOO 

t  I  I 

^  s 


UJ 

>  1- 

0  3 

0 

^  2 

z 

3  z 

cn 

?  < 

UJ 

5 

2  X 

K 

UJ 

H 

U- 

cn 

E  2 

0 

> 

cn 

0 

cn 

►— 

cn 

K  0 

H 

UJ 

<  0 

cn 

H 

X  0 

IT) 
Ul  2 

^0 

□c 

-1 

UJ 

H 

0 

< 

< 

X 

LU 

H 

< 

cn 

a 

X 

CO 

oc 

< 

X 

0 

s 

0 

00 

0 

_J 

3 

< 

1-  ** 

^  z 

X 

cn  < 

cn 

Q 

Z  X 

H- 

> 

0  K 

X 

cj  cn 

I  I 


(T 

LU 

CL 

O 

cc 

CL 


ZD 

O 

cr 

X 


< 

o 

!5 


3 

< 

2 

UJ 

< 


< 

Q 

I 


CE 

UJ 

H 

3 

CL 

S 

O 

o 


fO 


CL 

O 

a 


cn 

X 

UJ 

2 

o 


o 

z 

o 

o 


X 

3 

X 


Q 

LU 


CO 

UJ 


UJ 

X 

3 

H 

X 


X 

< 

X 

o 

o 

X 

CL 

z 

UJ 

s 

o 

UJ 

X 

cn 


a 

< 


niCROSTRAlN 


</) 

z: 

UJ 


o : 

LU 


UJ  * 


0) 

c 

0 

a 

T5 


:  LU 


Si 


E 

0 

0 

"D 

Hi 

c 

c 

0 

0 

U) 

0 

— 

— 

(/) 

c 

-t-) 

® 

0 

0 

c 

E 

c 

c 

0 

— 

— 

u 

E  ^ 

E 

•  — 

® 

0 

u 

0 

X 

— 

— 

Hi 

L 

Hi 

> 

~0 

“0 

— 

Hi 

D 

® 

L 

E 

L 

0 

0 

E 

0 

X 

® 

_ 

> 

— 

Ui 

D 

W 

:: 

■D 

a 

U 

c 

u 

0 

E 

L 

0 

L 

L 

0 

— 

c 

■  — 

X 

u 

CJ 

CJ 

i 

O 


* 

II 

(nm 

I 

I 


Ck:  M 


• 

1  uiuuBlS 

Z] 

ujuifr  gz 

• 

uiui  1  61 

uiuji  21 

uluj  gc  9 

P  uiui  0j  t 

IZ 

NIOHIA 

.  D  < 


I8S888888 
odw  ‘ssaais  3aniovai 


a: 

o 


o 

< 

Q- 


<0 

3 

I- 

< 

H* 

(O 


Si 

o 


■0 

0^ 


"0 
Hi  Hi 
W 

C  U 
0  Hi 
Q.X 
W  U 
Hi 

L  0) 

c 

0)  - 
0) 

0  -D 


“0 

c 

0 


c 

0 

—  c 

0 

0)  - 
-tj  y) 

c 

X  ® 


3  C 


w 

c 

® 

E 

0 

® 

a 

y> 

-D 

® 

o 

D 

a 

E 


-  W 


*D 

® 

(/) 

Hi 

W 

(A 

0 


C 

a  u) 

c 

T> 

0) 

0 

■— 

> 

0 

®  “0 

w 

c 

u 

u 

> 

— 

® 

® 

— 

-  *D 

■  — 

V 

0 

Hi 

E  C 

... 

05 

3 

9) 

u 

X 

0 

0  0 

•— 

c 

L 

"0 

X 

u 

c 

u 

Hi 

0  "0 

0 

0 

>“0 

0 

C  1? 

E 

c 

05 

•— 

— 

*D  ® 

0 

c 

0 

c 

L 

■  — 

D 

(A 

u 

Hi 

® 

L. 

u  C 

C 

L 

0 

— 

0  ® 

0 

0 

C 

0 

(A 

•  — 

0 

£ 

Hi 

0 

■  — 

(A 

0 

L 

c 

Hi 

c 

U) 

c 

c 

U 

V) 

0 

W  ~ 

0 

® 

— 

— 

•— 

“  a 

E 

■i 

0 

c 

E 

D 

0 

® 

(/) 

W  E 

■  — 

c 

■  — 

— 

D 

c 

L 

> 

0 

>  - 

L 

u 

E 

L 

~o 

® 

— 

L 

— 

0 

0 

E 

0 

® 

E 

<— 

0 

® 

0  “0 

a 

Cl 

0 

a 

(A 

E 

— 

c 

a 

c  c 

X 

E 

X 

® 

X 

0 

0 

< 

0 

<  0 

UJ 

M 

® 

Cli 

Hi 

Qi: 

u 

X 

• 

• 

■ 

• 

• 

• 

• 

